Kinetic Wave Turbulence 



Gregory L. Eyink 1,2 and Yi-Kang Shi 1 

(1 ) Department of Applied Mathematics & Statistics and (2) Department of Physics & Astronomy, 
The Johns Hopkins University, Baltimore, MD, USA 



Abstract 

We consider a general model of Hamiltonian wave systems with triple resonances, in the standard kinetic limit of 
a continuum of weakly interacting dispersive waves with random phases. We show in this limit that the leading- 
order, asymptotically valid dynamical equation for multimode amplitude distributions is not the well-known equation 
of Peierls (also, Brout & Prigogine and Zaslavskii & Sagdeev), but is instead a reduced equation containing only 
a subset of the terms in that equation. Our equations are consistent with the Peierls equation in that the additional 
terms in the latter vanish as inverse powers of volume in the large-box limit. The equations that we derive are 
the direct analogue of the Boltzmann hierarchy obtained from the BBGKY hierarchy in the low-density limit for 
gases. We show that the asymptotic multimode equations possess factorized solutions for factorized initial data, 
which correspond to preservation in time of the property of "random phases & amplitudes". The factors satisfy the 
equations for the 1-mode probability density functions (PDF's) previously derived by Choi et al. and Jakobsen & 
Newell. Analogous to the Klimontovich density in the kinetic theory of gases, we introduce the concepts of the 
"empirical spectrum" and the "empirical 1-mode PDF". We show that the factorization of the hierarchy equations 
implies that these quantities are self-averaging: they satisfy the wave-kinetic closure equations of the spectrum and 
1-mode PDF for almost any selection of phases and amplitudes from the initial ensemble. We show that both of these 
closure equations satisfy an //-theorem for an entropy defined by Boltzmann's prescription S - kg log W. We also 
characterize the general solutions of our multimode distribution equations, for initial conditions with random phases 
but with no statistical assumptions on the amplitudes. Analogous to a result of Spohn for the Boltzmann hierarchy, 
these are "super-statistical solutions" that correspond to ensembles of solutions of the wave-kinetic closure equations 
with random initial conditions or random forces. On the basis of our results, we discuss possible kinetic explanations 
of intermittency and non-Gaussian statistics in wave turbulence. In particular, we advance the explanation of a "super- 
turbulence" produced by stochastic or turbulent solutions of the wave kinetic equations themselves. 

Keywords: Keywords here. 



1. Introduction and Outline of Main Results 

The theory of wave turbulence has traditionally focused on the wavenumber spectrum, which is expected to satisfy 
the wave kinetic equation in the limit of a continuum of weakly interacting, phase-incoherent waves HI El- However, 
recent very interesting works of Choi et al. (2004,2005) ffl 16) and Jakobsen & Newell (2004) have studied 
higher-order fluctuations of wave amplitudes by deriving equations for probability density functions (PDF's) in this 
same kinetic limit. These works generalized the results of a pioneering study of Peierls (1929) [8], who obtained a set 
of evolution equations for multimode PDF's of energies of phonons (sound waves) in anharmonic crystals. Cf. also 
Brout & Prigogine (1956) [9]. Such equations were first discussed in the context of wave turbulence by Zaslavskii & 
Sagdeev (1967) [ 10 |, for some special 3-wave systems. The recent works cast new light on some outstanding problems 
of wave turbulence. On the one hand, the papers (3] |U [5] [6] |7j have made an important contribution to the founda- 
tions of wave kinetic theory, by showing that the "random phase" and "random phases & amplitudes" properties of 
initial wave fields are preserved in time by their multimode equations^ These developments are surveyed in a recent 



'Strictly speaking, such assumptions on initial data should be unnecessary for the validity of wave kinetics. RPA properties and, even more, 
Gaussian statistics should be dynamically established locally in physical space on the fast linear wave time-scale, due to the superposition of 
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monograph of Nazarenko lfl3l . Although there has been some recent progress lfT4l . we are still far from a rigorous 
mathematical derivation either of the wave kinetic equation or of any of its predictions. Another important result of 
G]|6ll3 is a closed equation for the 1-mode PDF's of the independent wave amplitudes. These equations can have 
as solutions 1-mode PDFs that are far from Gaussian, especially when the equations are supplemented with boundary 
conditions and additional terms to represent strongly nonlinear process, such as wave-breaking, that lie outside the va- 
lidity of weakly-interacting wave kinetics IUJQ3). These results show promise to explain observations of intermittency 
and anomalous scaling in some recent experiments lfT6l[T5l and simulations ifTTl of gravity-capillary wave turbulence, 
especially for wavenumber regimes where such intermittency was not previously expected lfl8l 

Unfortunately, the Peierls-Brout-Prigogine-Zaslavskii-Sagdeev (or simply, Peierls) equations for multimode PDF's 
are not the leading-order, asymptotically correct equations in the standard limit of a continuum of weakly interacting, 
incoherent waves. It was not observed in ||5] [71 that many of the terms in the Peierls equation vanish as inverse 
powers of volume in the infinite-box (continuous wavenumber) limit. Using the same methods as theirs but with a 
correct estimation of the terms, we shall obtain a novel set of multimode equations simpler than the Peierls equation 
On the other hand, we shall show using these new multimode equations that the 1-mode PDF equations obtained in 
[6l [7J are correct, under the assumption that initial wave modes have statistically-independent random amplitudes 
as well as phases. More generally, we shall classify all possible realizable solutions of the new multimode equations, 
for initial conditions with independent (and uniform) random phases but with no assumption on the statistics of initial 
wave amplitudes. We find that the most general solutions correspond to " super- ensembles" of solutions of the wave 
kinetic equations with random initial conditions for the spectrum or with random forcings. This represents another 
possible mechanism for intermittency and non-Gaussian distributions by a "super-turbulence" of chaotic or stochastic 
solutions of the kinetic equations. As we discuss in detail, existing results on linear stability of Kolmogorov cascade 
solutions of the wave kinetic equations do not rule out transition to such "super-turbulence." In fact, this possibility 
was anticipated in a review of the stability theory by Zakharov et al. (1992), (IT), section 4.2.2), who referred to 
this possibility as "secondary turbulence." The results of the present work show that this is the only possibility for 
explaining intermittency and anomalous scaling of wave turbulence strictly within the wave kinetics framework, fj 

The invalidity of the Peierls equation as a leading-order description will be shown below by specific, detailed 
calculations, but it can be understood on the basis of simple, general considerations. It will be useful to present such 
arguments here as a preliminary to the concrete calculations. This will also permit us to give an overview of the 
important new concepts, tools and results of our paper. We must begin with a brief summary of the main results of 
|7j|6j for a general Hamiltonian dynamics with 3 -wave interactions: 

• <r _ V" 1 Y^ 1 jcr,cr,,o-2 o-, 0-2 /(0-1 (d(ki)+<r2w(k 2 )-c™(k))fr /I s ! 
k / j / 1 Hjufe "k, "k 7 e Oo-,k l+ o- : k 2 ,o-k (l) 

0"1 ,0"2 ki ,k2 

where the wave triplet nonlinear interaction ^'^jf 3 = 3;cr 1 // k | T k ^' cr3 in terms of the coefficient of the cubic term in 
the Hamiltonian and where the summations range over wavenumbers k€ jZ'' and a degeneracy index cr = ±1, with 
L the side-length of a periodic box containing the wave system. (For a systematic discussion of the notations used in 
this paper, see section [TT| below.) The analysis of Q[6) is based on action-angle variables {J^, tp^) for the linear wave 
dynamics, defined by a k = ^[T^e" :r ' ek . Note that the standard spectral density is related to the action variables as 

n(k)= lim(^)V k > (2) 

in the infinite-volume limit, with total wave action per volume given by N = J d d k n(k). Here the average (■) is over 
the ensemble of initial amplitudes and phases. In the large-box limit (L — > oo) followed by the weak nonlinearity limit 
(e — > 0), the spectrum is expected to satisfy the wave kinetic equation: 

n(k,r) = 36tt ^ f d d k 2 d d k 3 \H^\ 2 5{a ■ w(k))<fV ■ k) 

Z=(— 1,0-2,0-3) 



dispersive waves trains arriving from distant, uncorrelated regions of space. See II lll2l or, more rigorously, 1121 . Such properties are assumed 
valid at the initial time only to simplify the derivations. 

2 We say "strictly" because the approach of [5 6| invokes strongly nonlinear effects at sufficiently large wave amplitudes, outside the kinetic 
description. Our remark applies only to leading order in the small parameters e and 1/L, whereas higher-order corrections could also in principle 
lead to (weakly) non-Gaussian statistics in either physical or Fourier space. 
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<{n(k2, r)«(k3, r) - cr 2 n(k, T)n(k^, t) - <X3n(k, r)n(k2, r)} 



(3) 



where t = e 2 f is the nonlinear time. See ifTl fTSl . To obtain information about the statistical distribution of amplitudes 
and phases beyond the spectrum ( |42| ), the authors of (5] |6) exploit a generating functional defined by 



Zl(A,h) = ( ex P 



Z iJ-kJk 



+ 'Hkfk 



V k 



(4) 



In the same kinetic limit (L — » oo, e — > 0), for an initial distribution with independent, uniform phases, this generating 
function is claimed to satisfy the Peierls equation in the form 



dZ 
dr 



(5) 



k,<r 



with <x = (cri, cr 2 , <r 3 ),k = (k;,k2,k3) and with 5t(w) an approximate delta function of width 1/T in frequency space 
for some time T intermediate between the wave period and the time t - e~ 2 r. From this equation, EE) derived results 
on the higher-order fluctuations and statistical distributions of the wave mode amplitudes. 

To see that the results, as stated above, cannot be asymptotically correct, consider the limiting behavior of the 
generating functional Q when the field a k is assumed to be "RPA", i.e. to have amplitudes and phases for distinct 
Fourier modes given by mutually independent random variables. In that case, the generating functional must factorize 
for all independent variables, in particular for the amplitudes, as: 



ZdAp. = 0) = (exp 



Z //Ik A ) = Y\ ( ex P ( //l k4) )> 

> k /' k 



(6) 



where each factor is a 1-mode generating function Z^i^k', k) = / exp (//IkA) )• However, the definition of the spec- 
trum implies that the rescaled variable 

J k := (L/2n) d J k (7) 

must be 0(1) in the limit as L — > oo. For the 1-mode generating function this gives, employing the cumulant expansion, 
the result: 



Z (1) (^ k ;k) = (exp 



//l k 



In 



d \ 
Jk 



) = ex p z 



V /2tt 



pd 



where (J^Y is the pth-order cumulant. Inserting back into the product d6|, one obtains 



ZdA,H = 0) = exp 



Now assume that /l k = /i(k) for some smooth function /i(k). Since for L — > oo 



2;r 



Z^ p (k)<^> e 



k A"(k){jjy, 



(8) 



(9) 



(10) 



the contributions of pth-order cumulants are 0(L d{p Only the p — 1 contribution «(k) = (7 k ) survives for very 
large L 

Zl(A,M = 0) = ex p(' ^ ^( k )"( k ) + 0(L _d )| L^oo'exp^J" A i(k)n(k)j . (11) 

Thus, the generating functional Zl(A,/u) of an RPA field is completely determined by its spectrum n(k) in the limit as 
L — > oo and it then contains no information about higher-order fluctuations. This result, which we have derived here 
assuming existence of all higher-order cumulants, can be proved under much weaker assumptions (see section 2). 
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There is a simple interpretation of the above limit in terms of the empirical spectrum, defined as 



— 2 A/(k-kO. (12) 

For an RPA field, this quantity is a sum of a large number of independent variables]^] It is exactly analogous to the 
"Klimontovich density" or empirical 1 -particle density in the kinetic theory of gases IT191 . Unlike the usual spectrum 
n(k), the empirical spectrum is a random variable that incorporates the information about amplitude fluctuations. The 
amplitude generating function ZlW - Zl(^ /-f = 0) is just the characteristic functional of the empirical spectrum: 

Z L {X) = (exp (i j d d k i(k)« L (k)J| . (13) 

The previous limiting result for Zl(A) is mathematically equivalent to the statement that the empirical spectrum has a 
deterministic limit for L — ♦ oo which is the usual spectrum: 

lim n £ (k) = «(k). (14) 

L— >oo 

This limit is a probabilistic law of large numbers for the empirical spectrum. There is a precisely analogous law of 
large numbers for the empirical 1 -particle distribution in the low-density limit for the kinetic theory of gases, as first 
shown in 1975 by O. Lanford 12011211 . 

The above results should hold not only for RPA fields at initial times, but also for the evolved fields in the kinetic 
regime of wave turbulence, if the RPA property is propagated in time as expected. For example, [5 1, section 5, and 
ifTSl . section 11.4, argue that the RPA property will be preserved in a somewhat weaker form. As we shall discuss 



in detail below (see eq.(56i), these weaker forms still suffice to derive the limiting exponential expression (Hi for 
the generating functional. However, the Peierls equation does not have such exponential solutions! Thus, the current 
claims in the literature are not consistent. We shall see that when carried out carefully, with due regard to the scaling 
of various terms with e and L, the standard analytical methods of wave turbulence yield not the Peierls equation in the 
limit as first L — > oo then e — > 0, but instead the following simpler equation: 



Z[A,p] = - 36^,0 Yj [ d dk i dd h d % <5V • k)<5(£ • w(k))|Hjf| 2 

cr=(-l,a- 2 ,<x 3 ) 



(15) 

.„ J 6 2 Z 6 2 Z 6 2 Z 

X /l(ki ) < <X 2 CT 3 

\6A(k 2 )6A(k 3 ) 6A(ki)6A(k 3 ) 6A(ki)SA(k 2 ) 

This equation will be justified in detail in this work. Here we note only that it indeed differs from the Peierls equation. 
It contains only a subset of the terms in the Peierls equation, lacking in particular all the terms involving third- 
order derivatives in A. The important property of the above corrected equation is that it has the expected solutions of 
exponential form 



Z[A,fi,T] = exp\i j d d k A(k)n(k,r)l (16) 

where n(k, r) solves the classical wave kinetic equation Q, if the initial conditions Z[A,p.,T = 0] are also of this 
exponential form. The property of preserving exponential solutions implies a law of large-numbers for the empirical 
spectrum n/,(k, r) at times t > and is equivalent to the "propagation of chaos" by the BBGKY hierarchy in the 
kinetic limit for gases. 



The reader may worry that all terms appear to vanish in this summation except the single term for ki = k! However, it is important to keep in 
mind that distributions (generalized functions) like the Dirac delta function are only meaningful when smeared with smooth test functions. In that 

case, J d d k tfi(kjii L (]i) = IL^^jd ■'ki^O 4 ') and this sum indeed contains a large number of terms in the limit L — > oo. More physically, the 
empirical spectrum must be "coarse-grained" by integration over cells A that partition the wavenumber space. This corresponds to replacing the 
smooth test function </>(k) by the characteristic function ^(k) which = 1 for k 6 A and = otherwise. 
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The analogy of kinetic wave turbulence with the kinetic theory of gases is in fact quite close. As we shall show 
below, the above equation for Z is equivalent to a hierarchy of equations for M-point correlation functions of the em- 
pirical spectrum, which is exactly analogous to the "Boltzmann hierarchy" obtained from the BBGKY hierarchy for 
low-density gases ll20ll2D . Just as the Boltzmann hierarchy has factorized solutions for factorized initial conditions, 
so does the kinetic wave hierarchy for all multi-point spectral correlation functions. The factors in both cases solve 
the relevant kinetic equation, which satisfies an //-theorem corresponding to positive entropy production. For both 
hierarchies the general solutions without assuming factorized initial data are "super-statistical solutions" which repre- 
sent statistical mixtures of factorized solutions. This was first discussed for the Boltzmann hierarchy by Spohn [22 J. 
For wave turbulence these "super-statistical solutions" of the spectral hierarchy correspond to ensembles of solutions 
«(k, t) of the classical wave kinetic equation with random initial conditions «o(k). We shall discuss in this work both 
the mathematical derivation and the physical relevance of the "super-statistical solutions". These are a possible source 
of non-Gaussian statistics and intermittency not widely appreciated in the wave turbulence literature. 

Previous attempts to calculate intermittency effects within wave kinetic theory have employed the Peierls equation 
or its reduced forms for the PDFs of a finite number of modes |5]|6]|23). If the Peierls equation is not asymptotically 
valid in the usual wave kinetic limit (first L — > oo, then e — > 0), then does this mean that the previously claimed 
results for the finite-mode PDF's are also not leading-order asymptotic limits? The answer is yes, but fortunately only 
partially. The joint PDF of a set of M wavenumber modes , Jk M is characterized by its Fourier transform, the 
M-mode generating function: 



( M 

r(M) 



Z w (A u ...,Am,H = 0;ki,...,k M ) = exp 



(17) 



This is the same object which was introduced in the monograph of Nazarenko [13], section 5.4, except that we 
consider only a fixed number M of the total number N = 0(L d ) of Fourier mode^] We find that these objects do 
not individually satisfy closed equations but instead satisfy a hierarchy of equations that link M-mode functions to 
(M + 1)- and (M + 2)-mode functions: 



M „ 

36/tt 2 Z J d ^d% S d {a ■ ±pS(a ■ co(k.))\H^\ 



M 

z (M) 

d g 2 Z ( M+ 2) dZ {M+l) \ dZ (M+l) \ ) 

{(Aj + A;—) _ _ L _ -0=2/1;— — - -^-a,-= — L . 

1 J dAj dA 2 dA 3 ^3=0 1 dA 3 8Ai l *=° dA 2 dA^-= 0> 



(18) 



with the notations <x = (-1, cr 2 , CF3) and k^. = (kj,k 2 ,k 3 ). The equations in this hierarchy can also be obtained from 
the Peierls equation, by setting all except M of the A's equal to and by retaining only a subset of terms. In particular, 
the terms in the Peierls equation which contain products A^A^ for k + k' are found to vanish in the limit L — > 00, 
when only a finite number M of k-modes have A^ + 0. The above hierarchy for the M-mode generating functions is 
equivalent to a hierarchy for the joint PDF's !P (M) (s 1 , sm', ki , k^) of M Fourier modes ki , k^: 

M „ 

pm + y = 0; (19) 



where, with the same notations as above, 



T) M) = - 36nsj J d d k 2 d d k 3 6 d (z • kj)6(g: ■ w(k .))|//|| 2 



tr=(-l,CT2,0~3) 

dp (M+2 > nm 

ds 2 dsi, — s M , s 2 , Si)s 2 sj, ^> 

dsj 



+ a 2 J ds 2 r ( - M+1 H Sl ,...,s M ,-S3)s i +^3 J d-s 2 ^ M+1 \ Sl ,...,s M ^ 2 Ts 2 ]. 



It is easy to see that if all N modes are considered, then Z^ N> diverges exponentially in the limit L — > 00 and. in particular, cannot converge to 
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a solution of the Peierls equation. Only the quantity {jf\ log-Z'^' can have a finite limit. We shall study this in detail in a following paper 1241 . 



Previously proposed closed equations for multi-mode equations, e.g. for M — 2 in [23], are not asymptotically exact 
in the usual kinetic limit for wave turbulence but require additional ad hoc closure assumptions. 

On the other hand, the one-mode equation that was derived in J6j for the PDF P(s, t; k) = !P (1) (s, t; k) of a 
single amplitude A : 



d_ 
Fs 



dP 
7,15 ~ds + 



(21) 



is obtained from the above hierarchy as a leading-order asymptotic equation, for factorized solutions of the form: 



P <M) ( Sl , sit, t; k, , k M ) = f] P(s m , r; k m ). 



(22) 



m=l 



The coefficients rj^, that appear also in the kinetic equation hk — % - Tk"k are obtained from the condition 
«k = J ds sP(s; k) and the standard wavenumber integrals over products of n. The one-mode equation is thus a 
nonlinear Fokker-Planck equation in the sense of McKean [25]. Factorized solutions of the form (22i are obtained 
for factorized initial data, corresponding to RPA fields. The hierarchy thus preserves in this sense the statistical 
independence of amplitudes. Just as for the spectral hierarchy, this factorization property of solutions implies a law 
of large numbers for empirical 1-mode PDFs 



1 



N L (A) 



2 S(s - A) 



(23) 



with A c R d a subset of continuous wavenumbers and with Nl(A) the number of elements in ~Z d n A. Also as 
for the spectral hierarchy, the general solutions of the PDF hierarchy are "super-statistical solutions" that correspond 
to ensembles of solutions of the 1-mode equation (21 1 with random initial conditions Po(s;k). As we shall discuss, 



such ensembles may be physically relevant to explain non-Gaussian statistics in weakly nonlinear wave systems if the 
solutions of (21 1 themselves become random, either through fluctuating driving forces or internal chaos/turbulence. 



The detailed contents of this work are as follows. The introduction concludes with a brief sec tion [TT| summ arizing 
our model assumptions and notations. The next section 2] discusses the probabilistic properties of "RP" and "RPA" 
fields. The principal results of the paper are in section [3 which derives the multimode equations and discusses their 
basic properties. This is done first for the spectral generating functional and correlation functions (section 3.1 1 and 
next for the PDF generating functions and correlation functions (section [3~2| . Many technical details are provided in 
|Appendix A| and | Appendix B| Finally, we discuss in section |4]possible kinetic theory mechanisms of intermittency 
in wave turbulence. The conclusion section [5] summarizes our results. 



1.1. Model and Notations 

We summarize here briefly the notations employed in our work. We consider a system consisting of a complex 
wavefield u(x,t) in a (/-dimensional periodic cube with side L. As in 17][6], we assume that there is a maximum 
wavenumber £ max , to avoid ultraviolet divergences. This can be achieved by a lattice regularization with spacing 
a = L/M, for some large integer M, so that £ max = n/a. The location variable x then ranges over the physical space 



A L = aZt, 



(24) 



with the usual notation Z M for the field of integers modulo M. This space has volume V = L d . The dual space of 
wavenumbers is 

Al = 2 4K (25) 



with fcmjn = 2tt/L. The total number of modes is N - M d , so that spatial volume V = Na d . We use the following index 
notation 

m(x) <t = +1 

M*(x) <T = -1 



K°"(X) 



(26) 
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for u and its complex-conjugate u* , following [7]. Likewise, we adopt their convention for (discrete) Fourier transform 

A ff (k) = - J] f) ex P(- /crk ' x ) ( 27 ) 

xeA t 

so that A + (k) and A~(k) are complex conjugates. Notice that this quantity converges to the continuous Fourier trans- 
form ji f[ L yt d d x u°~(x, t) exp(-/crk • x) in the limit a — > 0. The discrete inverse transform is 

u a (x) = A<J ( k ) exp(/crk ■ x). (28) 

kGA^ 

The dynamics is assumed to be canonical Hamiltonian with a cubic Hamiltonian density (energy per volume) 
describing 3-wave interactions: 

H = H +8H = Yj ^)\K\ 2 + e Yu H k A k,' A k" A k^-k,o- (29) 

k a,k 

As in [7 1, we denote triplets of variables by the notation <x = (cry , cr 2 , cr 3 ), k = (kj , k2, k-}) and define the dot product 
<x • k = crjki + <r 2 k2 + o"3k3. The Kronecker delta function in the triplet interaction term enforces the symmetry of the 
dynamics under the group of (discrete) space translations. The interaction coefficients are furthermore assumed to 
satisfy 

C3Q) 

for any permutation n e 53. The first condition guarantees the reality of the Hamiltonian and the second can always 
be assumed without loss of generality. The Hamiltonian equations of motion for the Fourier coefficients are 

k o- 1 ,o- 2 ki,k2 

where I-x^'Q = 3;cr// k ^'° k 1 ^ 2 . It is useful to introduce action-angle variables for the linear dynamics, = |A£"| 2 and 
(f>b = crarg(A^), so that A^ = VA*/^, where ifr^ = expiiipb). In these variables, the Liouville measure fi conserved by 
the Hamiltonian flow can be written variously as 

d/i — J~~[ dqkdpk - Y\ -dA^dA k = J~~[ dJ^dipk (32) 

k k ' k 

where the canonical momenta and coordinates are given by real and imaginary parts of A£ = -j^(Pk + icqt). 
As usual in wave kinetic theory, we introduce the "interaction representation" 

al = A ^ e -io-co(k)t (33) 

which removes the rapid wave oscillations. The Liouville measure is invariant under this transformation because of 
the rotation-invariance of the Haar measure on t/r e S . We shall often use the shorthand notation a\ = together 
with^i := Zcn.k, and 

w 23 : = -o"i^(ki) + cr 2 w(k 2 ) + cr 3 w(k 3 ), (34) 

^23 ' = ^o"ik], o- 2 k2+cT3k3 ■ (35) 

The dynamical equation of motion can then be written succinctly as 

a^e^Lm^^'^. (36) 

2,3 
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2. Fields with Random Phases and Amplitudes 



It is often assumed in derivations of wave kinetic equations that initial fields have Fourier coefficients with random 
(and statistically independent) phases and amplitudes. Furthermore, these properties are expected to be preserved in 
time, in some suitable sense, in the wave-kinetic limit. See the discussions in [7| and, particularly, ]5]|6). It is 
therefore important to review here the definitions and properties of such "RPA" field^] Our discussion shall be 
largely complementary to [6|, emphasizing the probabilistic characteristics of such fields. 

Consider then N complex-valued random variables a k , k e A* taken to be the Fourier coefficients of a random 
space field: 

«(x) = fl k exp(/k • x). (37) 

keA- L 

Here «k corresponds to = in the previous section (no distinction need be made between the two at time t = 0). 
It will be very important in what follows to work with normalized variables 

_ / L f 2 

° k \27r/ flk ^ 

which are assumed to remain finite in the large-box limit L — > oo. As we shall remind the reader below, this normal- 
ization is required so that the spectrum of the random field is well defined in that limit. It is convenient to write the 

complex variables in polar coordinates a k = VA<? Vk (action-angle variables) with the normalized action A = ( ^) Jk- 
We use the notations and ^ for possible values of the random variables £ K + and t^k = e ,iPk e S 1 . The Liouville 
measure of the previous section becomes 

dfi(s,0=Y]ds k 1 ^ (39) 

k£A* 71 

suitably normalized. We define the A^-mode joint probability density function P (N \s,^) with respect to Liouville 
measure, such that the average of the random variable f(J, if/) is given by 



= J dn{s,$<pW{s,&f{s,0 (40) 



where the integral is over (s, £) in the product space (M + ) A ' x (s . 

Following JT] |6] we define u(x) to be a random-phase field (RP) if tfr^ = e m for all k e A^ are independent and 
identically distributed (i.i.d.) random variables, uniformly distributed over the unit circle S 1 in the complex plane. In 
terms of the joint PDF, this is equivalent to the condition that 

P w (,s,£) =P m (g) (41 ) 

independent of It is easy to see that an RP u(x) is a homogeneous random field on A£, statistically invariant 
under space-translations by the finite group a% d u . This follows from the rotation invariance of the Haar measures 
\d^\/2n on the phase variables and by the representation of the translations u(x) — > u(x + r) for r € aZ d M as phase- 
rotations: i/'k — > iAk<?' k ' r . In the limit L — > oo the field ml(x) defined with appropriately chosen 7k,L will converge 
to a homogeneous random field m(x) invariant under translations by all 1 . The standard definition of the spectrum 
n(k) = lim^ 00 (L/2w) fl '(|flkxl 2 ) implies that one must choose 

lim(7 kt>L ) = «(k), (42) 

for k € A* = [-£ max , +k may ~\ d , where k L = ^(modM) • j- e A£ converges to k as L = aM — > oo (for fixed a). In this 
case, ul{x) converges in distribution as L — > oo to a homogeneous field u(x) with spectrum «(k). 



5 More accurate acronyms would be "IUP" for "independent uniform phases" rather than "RP" and "IPA" for "independent phases & amplitudes" 
rather than "RPA". However, we shall stick here to the abbreviations already employed in the literature. 
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Again following ||5]|6), we define u(x) to be a random-phase and amplitude field (RPA) if m(x) is RP and if also 
are mutually independent random variables for all k e A^. This is equivalent to the factorization of the iV-mode PDF 
into a product of 1-mode PDF's: 

^(i)=f]%;k). (43) 

k£A« 

All homogeneous Gaussian random fields are RPA. Conversely, for any sequence of RPA fields satisfying condition 
( 42 1 the spatial field «l(x) converges in distribution to the homogeneous Gaussian field with mean zero and spectrum 
n(k) as L — > oo. See [26 1 for related rigorous results. Here we note only that 

\d/2 

IJ KL exp(/k • x + i^k) (44) 



u L {x) 



(C2 

v ' k£At 



is a sum of N independent variables scaled by 1 / y/N. It is important to emphasize that the Fourier coefficients o^.l can 
remain far from Gaussian in this limit. In physical space also there are non-vanishing cumulants for large but finite L. 
A curious property of RPA fields for finite L is the slow decay of their cumulants in certain directions. For example, 
it is easy to calculate the 4th-order cumulant as 

(u(x l )u(x 2 )u t (x 3 )u*(x 4 )) - [<m(xiV(x 3 )Xm(x2)m*(x 4 )> + (1 <-> 2)] 

Iry \2d 

keAI 



|) j d d k [<^>-2</ k > 2 ] e ** 



+X2-X3-X4) 



(45) 



asymptotically as L — > 00. As expected, this goes to zero as 0(L~ d ) for L — > 00, but for finite L the cumulant is 
constant on the hyperplane Xi + X2 = X3 + X4, even when the pair of points Xi, X3 are separated very far from the pair 
X2,X4, for example. This is in contrast to some more traditional derivations of wave kinetic equations, which assume 
rapid spatial decay of higher-order cumulants ll27l . 

A most important result for RPA fields is that the empirical spectrum 

\d 

^ " (46) 

k,eA; 



(2 \^ 
T ) 2 4,L^(k-k!), k e A* 



converges under the condition ( |42| > to the deterministic spectrum n(k) with probability going to 1 in the limit L — > 00 
(weak law of large numbers). It is worth sketching the simple proof. Note for any continuous function A on A* that 



d d k A(kJn L (k) 



In 



2 ^(ki)-/i 



ki,L- 



(47) 



kieA; 



Crucially, 



In 



ki£A! 



£ AfkJnfkJ 
2 /l 2 (ki)<|/ kl -n(ki)| 2 > 

kieA^ 



^^ 2 (kX|7 k -«(k)| 2 > = 0(L- rf ) 



(48) 



under the modest assumption that J d d k n(k)| 2 ) < 00. Since also lim^oo i^f) 2a, eA^ ^(ki)«(ki) = J d d k A(k)n(k), 
the L 2 -convergence follows 

J d d k A(k)n L {k) - j d d k A(k)n(k) 
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lim 

L— >oo 



= 



(49) 



and thus J d d k A(k)ni(k) converges in probability to J d d k /l(k)«(k) for every bounded, continuous /l^This is suffi- 
cient to infer that the amplitude generating function defined in the introduction satisfies [^] 

lim Zl(X) = exp [i j d d k A(k)n(k)\ (53) 

with n(k) the deterministic spectrum. The law of large numbers (LLN) derived above means that for RPA fields 
the empirical spectrum «£(k) coincides with n(k) at large L for almost every realization of the random phases and 
amplitudes, not just after averaging over these variables. 

Notice that for the above result one does not actually need the full independence assumption in RPA, but only 
uncorrelated amplitudes: 

<44> = <4><4>. ki*k 2 . (54) 

An even weaker and more general condition can be stated in terms of the M-mode correlation functions of the empir- 
ical spectrum, defined as: 

N ( r M \k u ...,k M ) = <ni(ki) • ■ -n L (k M )). (55) 



Note that (42 1 implies that lim^oo A/£ (k) = n(k). A careful examination of the previous proof shows that, in order 
to obtain the LLN for the empirical spectrum, it suffices that 

lim [Aff (k,,k 2 ) - JV^kO/^fe)] = 0. (56) 

This condition is the analogue of the Stosszahlansatz invoked by Boltzmann to derive his kinetic equation. Under this 
assumption, all of the M-th order correlations that exist will factorize in the large-box limit: 

M 

lim N ( L M \k u ...,k M ) = Ff «(k m ). (57) 

m=\ 

These observations go back to Lanford in his derivation of the Boltzmann equation GU1I2D . The results that we shall 
present below suggest that properties (42 1 and (56 1 for the wave field at the initial time, together with the RP property, 
are sufficient for the validity of the wave kinetic equation and a LLN for the empirical spectrum at positive times. 

RPA fields whose Fourier amplitudes possess the full independence property satisfy an even stronger law of large 
numbers for the empirical 1-mode PDF, which is defined as 

(2 \ d 
2 *(s-4)^(k-ki). (58) 
' k,£At 



technically, this implies convergence in the weak-* topology for «/.(k), n(k) as bounded, positive measures on A*. More physically, the 
function /l(k) may be taken to be a smooth kernel Gs(k - ko) of width r5 centered around a particular wavenumber ko. The coarse-grained spectrum 
«i,<s(ko) = f d d k G$(k - ko)«z.(k) can then be interpreted as the result of measuring the spectrum at ko with a finite resolution 6 in wavenumber. 
The weak-* topology implies the convergence of such coarse-grained spectra in the limit L — > oo, when arbitrarily many wavenumbers lie within 
distance 6 of ko. 

7 This is another standard result in probability theory, but, for completeness, we here recall the proof. Let Xi = J d d k /l(k)7i/,(k) and x = 
J A , d d k i(k)n(k). For e > 0, write 

\(e iXL - e ix )\ < (\e i( - XL ~ x) - 1| : \X L - x\ < e) + ^ Xl ~ x) - 1| : \X L - x\ > e) (50) 

For any S > 0, one can choose e so that \e' z — 1| < d when \z\ < £■ The first term is thus bounded by 6. By Chebyshev inequality the second term is 
bounded as 

(\ e m-x) _ n : \ Xl _ x \ 2P(\X L - x\ > e) < ^{\X L - x\ 2 ) -» (51) 

for fixed e as L — > oo. Thus, 

limsup|(e' Xi - e fa )| < S. (52) 

L— >oa 

Since 6 is arbitrary, it follows that lim^-^ (e lXl ^ = e' x . 
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Assume that the limiting random variables = lim/,^ M Jy Ly L of an RPA field exist and have PDF's P(s; k) which are 
continuous in k. Then an analogue of the previous argument implies that the random functions Pl(s; k) converge to 
P(s; k) with probability approaching 1 as L — > oo. This implies the previous LLN for the spectrum, since «z,(k) = 

poo ' — ' poo 

J ds sPi(s;k) and n(k) = J Q dssP(s;k). Although the "empirical PDF" defined in (58 i is mathematically very 
convenient, it is not a PDF for finite L. It is therefore more intuitive to use an alternative definition 



P L (s;A) = -±— V 6(s-J k ), 
N L (A) ^ 



(59) 



keAt nA 



for any open set A c A* and with Ni(A) the number of elements in A* n A. This quantity is nearly the same as 
j^j f A d d k Pl(s; k) for large L but it has the advantage that it defines a probability measure in s for each fixed A and L. 
Definition ( 59 1 also has a simple intuitive meaning, since it represents the instantaneous distribution of amplitudes of 
the large number of Fourier modes that reside in the set A for large box-size L. Under the same assumptions as above, 
it follows with probability going to 1 that 



lim P L (s 



; A) = — f d d . 
|A|J A 



k P(s; k) = P(s; A). 



(60) 



As before, strict independence is not necessary for the above results to hold. The essential property that is required 
is a factorization property of multimode PDF's, defined for kj , e A* by 



^ M) (s!,..., s M ;k!,...,k M ) = (6(si - Jk hLjj ) ■ ■ ■ 6(s M - A^))- 
The factorization property of the 2-mode PDF's for all pairs of distinct ki, k2 e A* 

lim [<pf\s , , S2 ; ki, k 2 ) - pV\ Sl ; k^fe; k 2 )] = 

suffices^] to derive the LLN for the empirical PDF and also the factorization of the multi-mode PDF's 

M 

lim !P^ M) (ii,..., s M ;ki,...,k M ) = M P(s m ;k m ) 



(61) 



(62) 



(67) 



for all integers M > 2 and distinct ki,...,kM e A*. The "asymptotic independence" property (62i is considerably 
weaker than strict RPA, permitting statistical dependence between Fourier modes at finite L. We shall discuss natural 



"microcanonical measures" in section 3.2.2 which satisfy the condition d62b but not the more stringent RPA condition 



In the following section we shall show that the above "generalized RPA" properties (57 1, (62 1 of initial conditions 
are preserved in time by the limiting kinetic hierarchies of wave turbulence. 



<?k (p(s,kXP L (s;k))\ ) = (64) 



8 Since the proof is so similar to the one given previously, we give just a few details here. Because of our assumptions on = limL-,^ Jk Li L it 
is easy to check that the average of the empirical PDF converges as L — > oo : 

x E ^ )(s; k ' )sd(k ' - k) ^ P(s ' k) - (63) 

Therefore, in order to prove that Pt(s; k) converges in probability to P(s; k), with convergence in the weak-* topology on regular Borel measures, 
it is enough to show that 

Urn i\J ds J c^k<p(s,\i)P L (s;k)- J ds J 
for any continuous function <f(s, k) on R + X A* which vanishes as s -t oo. Now a direct calculation of the above average shows that it equals 

The bracketed expression for k'j = k^ is found to be 

pf\s u s 2 ;K^0 ~ ifUnK^L^K) = S(s, - S2)P™Ui;k\) - P ( l\ Sl ;k' 1 )lf( S2 ;k' l ) (66) 

so that its contribution to the average vanishes as ~ (je) f A , d d k ^|v(/k, k) - {(p{Jk, k))|~^ for L — ► oo. Thus, to prove convergence, it must only 
be shown that the bracket term for k\ # k' 2 vanishes as L — > oo. This is obviously true for RPA fields, when the bracket term for k'j # kj is zero! 
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3. Multi-Mode Hierarchy Equations 



In this section we present our formal derivations of the multimode equations of wave kinetics and a discussion of 
their main properties. Our derivations are quite similar to those of Choi et al. (6), but our model is somewhat more 
general in appearance. The major difference between our analysis and theirs, however, lies in a correct accounting of 
the size of the various terms in the limits as L — > oo and e — > 0. 

We shall focus on the multimode generating function for finite box-size L, defined as in [6| by 

Zl&,H,T] = Lxv J]A k J k (T) Y]ifC(T)), (68) 



where A^ 6 K and jUk E Z for all k e A R . The time T is free for the moment but will later be chosen to be a time inter- 
mediate between the wave-period and the nonlinear time-scale. This generating function is calculated perturbatively 
in the nonlinearity parameter e at finite L, by expanding the solution of the dynamical equation 



:J]L 12 3fl2fl3^'^3 (69) 

into a power series 



2,3 



A straightforward calculation gives 



ai{T) = fli(0) + ea ( }\T) + e 2 af } (T) + 0(e 3 ). (70) 



af = a,(0), (71) 

a\ 1} = J] ^i23«f af A r (^ 3 )4, (72) 

2,3 

af = J] L^L 1Ai afaf af E T {u\ A5 , 0^)6^8% +(2^3) (73) 



2345 

We employ here the standard definitions l27l 

-*T r-T 



A T (x) 



I exp(ixt)dt, Ej{x,y) — I A,(x - y)exp(iyt)dt. (74) 
Jo Jo 



The terms in this perturbative solution of the equations of motion can be represented by a version of the Wyld diagram 
expansion |28 1. In this technique the various contributions are represented by tree diagrams, as illustrated in Figs. 1-3 
for the zeroth-, first- and second-order terms. In our conventions, similar to those of (6), a solid line labeled with an 
integer j represents a factor whereas a dashed line indicates the absence of such a factor. An arrow is added to 
a solid line to indicate crj = +1 ("source") when the arrow is pointed away from j and crj = -1 ("sink") when the 
arrow is pointed toward j. The triple vertex labelled 1, 2, 3 represents a factor L^'^Qe^'S^ with = +1 when the 
arrow points out of the vertex and cr\ — -1 when the arrow points into the vertex. The times at each vertex are ordered 
causally, with the latest times at the "root" of the tree, here labelled by 1 . When integrations are performed over all 
times from to T consistent with this ordering, then the various contributions to the perturbative solution result. 



Figure 1 : Zeroth-order terms and a? 



The generating function is obtained perturbatively by substituting ( 70 1 and expanding to obtain 

ZiU,jU, T] = X L {A,fi, T) + X* L {A, -n, T) 
12 



(75) 




Figure 2: First-order terras and a' 




Figure 3: Second-order term Oj 

with 

^R/i^} = ^R/z.O} + ( e AtJ *Wi + e 2 {J 2 +j2+jA+J5)]) J + 0(e i ), (76) 

keA* 

where, as in O, 

^Hn^E(* + A)^>, (77) 

k 1 ZJ l 

z k i 4 a ^ 

k 1 LJ \ 

* - < n I + ;>? - '> + jSJoWV w> 

k 1 LJ \ 

k 1*2 

77(0) ) 7(0) ( 2 a 2 ~ a 2 a 2 > a i a i )L- (81) 

The various contributions before averaging over phases can represented by diagrams, combining the tree diagrams for 
each of the factors by joining the trees with the same "root" indices. Each of the integer labels indicates an index to 
be summed over independently (except for the constraints imposed by delta-functions at the vertices). We illustrate 
this representation in Figs. 4-5 below for the first two contributions: 

The only contributions which survive the average over phases must have all phases summing to zero before av- 
eraging. This means that every a (0) factor must either pair with another factor a (0) so that their phases sum to zero 
or belong to a set of a (0) 's that pair with a ^°^ k factor so that the sum of all their phases is zero. The first we call 
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(1) JO)* 



Figure 4: First-order term 2 \(Al + jj~) a \ a 





Figure 5: Second-order term 2i + X\J\ ~ 3-7- )l<a I 



(11,2 



an "internal coupling", represented graphically by a solid line connecting the paired indices i, j which contributes a 
factor (5 -,+o- J ,o£k,,k J after phase averaging. The second we call an "external coupling", represented by joining all of 
the solid lines for indices i\, 12, i p to a blob • labelled a that represents the phase i/' k ° ,A ' k ° which contributes a factor 
(Jo-j + ... a - ip+Ala o Y\ P =\ "Skj.k,, a ft er phase averaging. This graphical representation is essentially the same as that employed 
in di). We illustrate the representation in Figs. 6 and 7 below for the terms in JTi and STi which turn out to give the 
leading -order contributions to those quantities as L — > 00. 



• 4 





Figure 6: Contributions to J~\ . 



Note that solid lines connected to external blobs have their wavenumbers "pinned" at the wavenumbers of those blobs, 
so that those wavenumbers are no longer summed over. For simplicity we often omit the labels of the blobs, since 
those play no important role other than the "pinning" described above. 

3.1. Spectral Hierarchy 

We now consider the first of the possible limits of as L — > 00 involving all of the modes. This leads 

to a set of equations for the spectral generating function and the spectral correlation functions. We first sketch the 



derivation of these equations, with more details in Appendix A and then analyze their basic properties. 
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Figure 7: Leading contribution to Jjj- 



3.1.1. Derivation 

The crucial observation which leads to our results differing from @ is that one must keep = 0(1), not 
0(1), in order to have a finite spectrum in the limit L — > oo. Thus, we take in the generating function of Choi et al. 

2n\ d 



A = | — | A, ^ = M(k), (82) 



where /l(k) is a smooth test function and, as before, ^ are integers. This leads to 



Zdl/i] = (exp i J] (j) A( V J * ft ^) ' (83) 

' keA; ^ ' J k£A* ' 



We consider the large-L asymptotics of the various terms in the perturbation expansion of this quantity, employing the 
following standard substitutions: 

From now on, we suppress superscript (0) when there is no confusion. 

Calculation of Si-' This quantity is represented by the graphs in Fig. 6, or analytically: 

J\ = ^ £123(^1 + 2j~) ,1^2+^2. 0^(^23)^23 



1*2^3 



+ J] ^122(^1 + yj-) V^1^2^i,l ( ^2+2cr2,oA(W22)^22 (85) 



1*2 1 



Note that o~\ — +1 in coefficient L123. Taking into account the wavenumber delta functions, there are two summations 
in the first term and one in the second. However, as noted by Choi et al. (2005), these sums contain only a couple 
of non-zero terms, and then only for special choices of the yu's. In their terminology, the sums are "pinned" by these 
choices of fi. Making the substitutions J, = (2n / 'Lf ^ 2 j ) and A\ = iA(ki) into the above gives the leading contribution 



(27r\^^ 1 \j J 

(2n\ cl ^ 1 1 J 

y 7! tm^i -y-^.i^^A^^ - L ~" 12 (86) 
L > m 2 V^T 

As we shall see, this term first-order in e gives a subleading correction in the limit L large but finite, larger than many 
of the terms that Choi et al. retained in their evolution formula for L — > 00. 

Calculation of Averaging over phases, the leading terms are contained in the contribution from Fig. 7 and a 
similar contribution with 2 and 3 interchanged. Because there are no external couplings, all fi's must vanish. Because 
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of the internal couplings, there are no sums over the wavenumbers 4,5 which are the same as wavenumbers 2,3. The 
final contribution to the generating function is: 

{e^ AkJk J 2 ) J = Stf J] ((A, + A\j l )j 2 J i e^ kAkJk ) ] |L£££ | 2 |A((r 2 w 2 + cr 3 aj 3 " wi^.^k, 



9Vo g(, 1 + , ; ^)_g_ x| ^ h4( 



cr 2 w 2 + cr 3 W3 - £k,, 0^2+0^3 



(87) 



Papers HID retain both of the terms in the first factor in (78 1, i.e. both A\ and A\J\. However, with the proper scaling, 
A\ = 0(1) while A 2 X J\ = A\Ji = 0(L~ d ). Thus, the second term is even smaller than the contribution from JTi 
in the limit as L — > 00 and should be neglected. This can be seen also substituting /, = (^f) d ^ 2 Ji and Ai = iA(k{) and 
taking the limit L — > 00 using ( 84 1: 



{e LkAkJk J 2 ) ] ~ -9id M , f d%d d k 2 d d h A(ki) + (^] i(ki) 



&l(kl) J 



<5 2 Z 



<w(k 2 )<M(k 3 ) 



(88) 



cr={-l,a*2,o"3) * 

The only surviving term in the limit as L — > 00 is the one proportional to /l(kj). Verifying our initial estimation, the 
term proportional to /l(k 1 )<5/5/i(k 1 ) is 0(L~ d ) and vanishes in the limit. 

A similar analysis may be carried through for the remaining contributions from ^3,^4 andJV Detailed discussion 
of all the terms is given in Appendix A where it is shown that only ^J 2 ,^J 3 ,^Js give 0(1) contributions in the large-box 
limit L — > 00: 



X(T) - X(0) 



+ £- 



f d d k x d d k 2 d d k 3 5 d (a ■ k)|A r (o_ • w(k))| 2 

2=(-l,CT2,0- 3 ) 

X/l(ki)flr — — — - 
1 5/t(k 2 )W(k 3 ) 

18/^o X cr 2 f d d k x d d k 2 d d h5 d (a-K)E T (0,oi- w(k)) 

' <r=(-l^,<r 3 ) 

xA(ki)|flff 



5 2 Z 



W(ki)^(k 3 ) 

1 k^ki-k! ^=(-1,0-2,0-3) ^ 



ErMki) + w(-ki).2 ' «(k))^-kXk 3 \ eZk ' kyk 

V„i [1 ^° Z f d d k 2 d d k 3 6 d (g:-k) 

1 k#±k, 



y'-A) +c<-> 3) 



2:=(-,o*2,tr 3 ) 

[A r (ff ■ w(k))A r (-<r • w(k'))] H^HXX* 



vax 



(89) 



The first term is that already found for J2, the second bracketed term is from J3, and the final term is from J5. Note 
that the expressions proportional to i5 /Jl ,i5 /J _ 1 , were missed in Q|6) (although they are larger than many extra terms that 
those authors retained in their final equations which actually vanish as L — > 00!) The terms proportional to S^id^, , 
indeed do not appear in the final equations, not because of the large-L limit but because they are nonsecular. 

We therefore consider the limit of weak nonlinearity. The limit is achieved by choosing the time variable T in 
the expansion (70 1 to lie between the wave period and nonlinear timescale 0(e~ 2 ) and by then taking e — > 0. For this 
purpose we use the following standard asymptotic relations for T — > 00 |27l : 

A T (x) ~ A(jc) = n8(x) + iP ( ~ J , E T (x; y) ~ A r (*)A r (y) ~ A(jc)A(y), 
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\A T (x)\ 2 ~2ttTS(x) + 2P^^-, E T {x;0)~A(x)(r-i-?-\, (90) 
The terms multiplied by S^o contain secular contributions proportional to T, while the terms proportional to S^iS^ , 



are nonsecular. We now use (75 1 to calculate Z(T) - Z(0), with the observation that it is changed by our replacement 
/ik = z'/l(k) into 

Z[A,», T] = X{A,p., T) + X*{-A, -/j, T). (91) 
Finally, replacing (Z(T) - Z(0))/T by Z and using time variable r = e 2 ?, one obtains 



d 
(It 



Z[A,fi,T] = -36inS^ Q ^ j d%d d k 2 d d k 3 5 d (a-k)S(a-co(k))A.Qs. 1 )\H~i 



£ 2 Z 6 l Z 5 l Z 



°~2 rfb ; ~°"3 rr^ ,7^, ; ^ (92) 



v cM(k 2 )<W(k 3 ) W(ki)M(k 3 ) W(ki)«(k 2 ) 
Equation d92j is the main result of this section. 



The validity of ( 92 1 does not, of course, require e infinitesimally small or L infinitely large, but just an e sufficiently 
small and an L sufficiently large (depending upon e). However, it is worthwhile to stress the precise conditions, largely 
following the prior discussions of [ 18 1 and 1 13 1, Chapter 1 1. In the first place, e must be so small that there is a large 
separation between the wave period and the nonlinear time 

1 , n(k) 

« e' 2 , „\, , , (93) 



a>{k) dn(k)/dT 



with the latter of order e y^. This is required in order to be able to find an intermediate time T so that ( 90 » and 
(Z(T) — Z(0))/T = Z are both well-satisfied, e.g. taking T = 0(e~ l ). As emphasized in [ 18 1, condition (93 i is almost 



never uniformly valid for all wavenumbers k, but is typically violated for either low or high k. This means that A(k) 
in |92| must be restricted to be zero for k outside the interval where |93| holds and it is furthermore assumed that the 
wavenumber integrations are sufficient local that no wavenumbers outside that range give a substantial contribution. 
The size of L is determined by the requirement that wavenumber summations over [-n/a, n/a] d can be approximated 
as continuous integrals. At the very least, it must be true that L » a. The most stringent condition seems to arise from 



the requirement that the approximate delta functions of width l/T in (90 1 contain a large number of frequencies. If 



Alo is the spacing of discrete frequencies, then one must have Aw •« l/T ~ 0(e). Estimating Aco = |Ak| • |VkO»(k)| ~ 
2 ^|Vk<y(k)|, the essential requirement is that L/|Vt<j(k)| » T, i.e. the time required for a wavepacket traveling at the 
group velocity to cross the box must be much larger than the time T or, conservatively, the nonlinear interaction time. 



This requires extremely large boxes in practice, unless e is only moderately small and (93 i only marginally satisfied 



3.1.2. Properties 

In this section we shall discuss the most basic properties of equation ( 92 1 derived in the previous section. As pointed 
out in [6 |, the factor <5^,o implies that the RP property of the initial conditions is preserved in time. Therefore, without 
loss of generality, we need only consider the characteristic functional for amplitudes, or Z[A, t] = Z[A,fi = 0, r]. Its 
evolution equation is 



Z[A,t] = -2,bin > I d a kid d k z d%6 d (<r ■ k)5(cr ■ co(k))A(ki)\H^ 2 



tr-(-l,CT2,o- 3 ) 

6 2 Z 5 2 Z 5 2 Z 



x{ — o- 2 — cr 3 — } (94) 

' <M(k 2 )<M(k 3 ) 6A(ki)6A(k 3 ) i 6A(ki)6A(k 2 ) j 



Hereafter we consider only this amplitude characteristic functional. Eq.( 94 1 implies a hierarchy of evolution equations 



for the M-mode spectral correlation functions defined in section|2] in the wave kinetic limit: 

N <M) (k u ...,k M ,T) = lim lim A/f?(ki,...,k 4ft e- 2 T). (95) 
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The hierarchy is easiest to derive by using the relation ( 13 i between Z[A, t] and the empirical spectrum, which implies 
that 

5 m Z[A,t] 



&Kki) • ■ ■ 6A(k M ) 



(96) 



,1=0 



By taking M functional derivatives of ( |94| and setting A = 0, one derives for each integer M — 1, 2, 3, ... the following 
equation: 



7V (M) (ki , k M , t) = 36n £ Yj f d d k 2 d d h5 d (v ■ kpd{a ■ co(k y ))|flg| a 

7=1 a;=(-l,er 2 ,a-3) 

[yV (M+I) (k, , k y _,, k i+1 , .... k„, k 2 , k 3 , t) - a 2 N (M+1) (ku k M , k 3 , r) - «7- 3 ^<" +1) (k 1 ,...,kj # ,k 2 ,T)]. 

(97) 

which couples the Mth-order correlation functions to the (M+ l)st. We shall refer to the above collection of equations 
for all M - 1,2,3... as the spectral hierarchy of kinetic wave turbulence. It is exactly analogous to the "Boltzmann 
hierarchy" derived by Lanford from the BBGKY hierarchy in the low-density limit |20 21 1. If the spectral correlation 
functions satisfy bounds on their growth for large orders M that allow them to uniquely characterize the distribution 



of the empirical spectrum, then the spectral hierarchy (97 1 is not only a consequence of the equation (94 1 but is in fact 
equivalent to that equation. 

d97l 



An extremely important property of the equations ( 94 1 or 
ticular, if the initial functional Z[A, 0] is of exponential form (p3b 



is that they possess certain exact solutions. In par- 
as follows for an initial RP field with uncorrelated 



amplitudes, then an exact solution of ( 94 1 is 



Z[A, t] = exp 



(i J d d kA(k)n(k, t)), 



(98) 



where n(k, r) satisfies the standard wave kinetic equation Q with initial condition n(k, 0) = n(k). This may be checked 
by direct substitution of (98 1 into (|94|l. Equivalently, with factorized Mth-order correlation functions (57 i as initial 



data, there is a solution of the spectral hierarchy equations (97 i also of factorized form: 



yV (M) (k 1 ,...,k M ,r) = f| n(k m ,T). 



(99) 



Note that such factorized solutions have only power-law growth for large orders M, so that they uniquely characterize 
the exponential characteristic functional (|98|. If it can be proved that solutions of the dynamical equations (|94| or (|97]> 



are unique for classes of initial data that include the forms (53 1 and (57 1, then the equations we have derived imply 
that spectral correlation functions initially factorized will remain so for t > 0. In this sense, therefore, the property of 
uncorrelated wave amplitudes is preserved in time by our equations. This is an exact analogue of the "propagation of 
chaos" property for the Boltzmann hierarchy, which implies that the Stosszahlansatz is propagated in time I201 I2T1 . 
The uniqueness of solutions of the linear functional equations ( |94] > and ( |97| ) is an interesting mathematical problem 
which we shall address in future work. 

The results above have an important implication. As follows from our discussion in section|2] the conditions (98 i 
or ( 99 1 imply a law of large numbers for the empirical spectrum at positive times. That is, with probability going to 1 
in the kinetic limit (first L — > oo, then e — > 0), it follows that 



n L {k, e~ 2 T) =* n(k, t), t>0 



(100) 



where n(k, t) is the solution of the wave kinetic equation ([3]). The interesting implication for laboratory and numerical 
experiments is that the wave kinetic equations will be valid for typical initial amplitudes and phases chosen from an 
RPA ensemble and not just for the spectrum averaged over the RPA ensemble. That is, the empirical spectrum is "self- 
averaging." This is the exact analogue of the law of large numbers derived by Lanford for the empirical 1 -particle 
distribution (Klimontovich density) in the low density limit for gases ll20ll2Tll . 
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In order to emphasize the close formal analogy of wave kinetic theory with the kinetic theory of gases, it is 
worthwhile to make here a few remarks about the role of entropy in both. It is well-known that the wave kinetic 
equation Q satisfies an "H-theorem" for the entropy defined, with Boltzmann's constant k B , by 



S[n] = k B d d k\nn(k). 



(101) 
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That is, iiS /dr > for general solutions of <j3j and dS jdr - for the thermal equilibrium solutions n eq (k) — k B T/a>(k) 

for fixed energy E[n] = J At d d kco(k)n(k). See Q, 
ly derived by Boltzmann for his kinetic equation, 
also follows from Boltzmann's prescription that 



at absolute temperature T, which maximize the entropy 

section 2.2.2. This is the exact analogue of the //-theorem originally derived by Boltzmann for his kinetic equation. 
As pointed out more recently by Spohn [29 1, the entropy 
S — k B logW [30 1, where W is the Liouville measure of the set of microstates (J, if/) consistent with the given 
"macrostate" defined by the prescribed spectrum «(k). More precisely, let A,, i = 1, .., P be a partition of A* and let 
Tn.A.^.L be the set of microstates {(J^, fa), k € A* } such that \n L (A,) - n(A,)| < rj, for i = 1, P where 



«l(A,) 



|A,| Ja, 



kn L (k), n(Aj) 



— f d d kn(k). 
|A/| Ja, 



TherQ 



(2 \ d P 
lim lim — In |r„, A ,,, L | = V |A,[ (In n(A ; ) + 1) . 

rj—*0 L—>oo \ L J *• — * 
> t i=\ 



(102) 



(103) 



The result for each individual cell is the same as the microcanonical entropy s(e) of a system of non-interacting 
harmonic oscillators or alternatively as the entropy of an ideal gas in a periodic box, thinking of Jk as kinetic energies 
and (fib as positions of the particles. Furthermore, defining |A| = max, |A,|, the partition may be refined by taking 
|A| -> 0, so that 

(2n\' 



lim lim lim 

|A|->0i7->0L->oo 



In |T„, A 



f d d 

J A' 



k (lnn(k) + 1). 



(104) 



The result agrees, up to constants, with ( 101 1. The above argument introduces the "microcanonical measure" obtained 
by restricting Liouville measure to the set T^a^l and normalizing by |r n> A,^,z,| to yield a probability measure. By 
adapting the arguments of ETI it follows that this "microcanonical measure" is a natural example which satisfies 



asymptotic factorization (57 1 but not the RPA property. 

The previous results may appear somewhat disappointing. Equations for the generating function Z,[A, r] like our 
( 94 1 have been proposed mainly in the hope of developing theories of intermittency of wave turbulence and of higher- 
order statistics of the wave-amplitudes. It might be concluded from the discussion above that the correct equation 
for Z[A, t] — i.e. our equation (94i — is equivalent to the wave kinetic equation <j3j and has no more physical content. 
However, this is not correct. To clarify this point, we now classify all of the realizable solutions of d94]i. That is, we 



characterize the solutions of the equation (94i for the most general possible initial conditions that can be physically 
attained in the kinetic limit, assuming RP but no form of RPA. For this purpose, note that the characteristic functional 
Z^zX^O) defined in ( 13 i is a positive-definite functional, i.e. 



2 CtcyZM- Aj> o)>o 



(105) 



for any set of n fields A\, ...,A„ and complex numbers c\, ...,c„. It is also normalized so that Z.l{A = 0,0) = 1. Since 
these properties are preserved under pointwise limits, the physical initial condition 



ZU,0)= lim Z L Wi) 

L— iOO 



(106) 



9 Note that r„ a, ;; .l is a Cartesian product set of the form ®? =i {(Jk, iftk), k e n A; : [?Tz,(A;) - n(A,)| < r/} and for each cell A, of the partition, 
«i(Aj) sb 2nI(A ) 2k£A,nA* (?k + <?k)' ^ e stated result then follows using the formula jpj^y R for the volume of a ball in dimension D = 2N(Af) 
of radius R = (nD) 1 ^ 2 , dividing by (2n) D ^ 2 [eq. |39| ], and applying Stirling's approximation for the factorial (D/2)\ as D — > oo. 
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for our eq.(94i also satisfies them. We must assume that the above limit exists, if eq.(94i is to have any validity at 
all, and we shall assume furthermore that the limiting Z.(A, 0) is a continuous functional of the A field^] By the 
Bochner-Minlos theorem, it therefore has the form 



Z(A,0) = j dp(n )exp(i j d d k A(k)n Q (k)) 



(107) 



for some probability measure p. E.g. see 1311 . That is, the initial condition Z.(A, 0) is a statistical superposition of 



exponential initial conditions of the form ( 53 i that arise from RPA fields. Since equation ( 94 1 is linear in 2, a solution 
for such superposed initial data is 



Z(A, t) = J dp(n ) exp (i j d d k A(k)n(k, t)) 



(108) 



where n(k, r) solves the kinetic equation with initial condition «o(k). This is the only solution if the uniqueness of 
solutions holds for (94 1. The conclusion of this argument is that the most general, statistically realizable solution^] 
of eq. ^4\ correspond to ensembles of solutions of the wave kinetic equation Q with random initial conditions «o(k). 

The above argument is a formal analogue of a rigorous result of Spohn for the Boltzmann hierarchy in the kinetic 
theory of gases Il22l . Under appropriate technical assumptions, his argument can be carried over to our spectral 



r M 
N {M) (k u ...,k M ,T) = J dp(n )Y]n(k m ,T), 



hierarchy (97 i, with the conclusion that its general realizable solutions are statistical superpositions of factorized 
solutions, that is, 

M 

(109) 

m— 1 

where again n(k, t) solves the kinetic equation with initial condition «o(k). We shall refer to such solutions of the 
spectral hierarchy ( 97 1 or of the equivalent equation ( 94 1 as "super-statistical solutions", since they correspond to 
random ensembles of solutions of the spectral closure equation. As we shall discuss at length in section |4] such 
"super-statistical solutions" offer a possibility to explain intermittency and non-Gaussian statistics previously little 
discussed in the wave turbulence literature. However, we shall first consider the alternative approach based on closed 
equations for the PDF's of the wave amplitudes. 

3.2. PDF Hierarchy 

We now consider a second possible limit involving only a fixed number M of modes a^, m = 1, M as the total 
number N — > oo. As before, one must keep — 0(1) for all modes. We thus define the joint characteristic function: 



Z ( L 1) (A u ...,A M ,p u ...p M ,T;k u ...,k M ) = (exp 



( M 



M 



/ m=l 



This is the generating function (68 1 of [6| with 



Ak "' = '(^r) A " u /k » = (^- /k - 111=1 A/ 



(HO) 



(HI) 



and all other A^ - 0. It also corresponds to the generating functional of [13|, eq. (5.15), with essentially the same 
definitions of A„, and 7™ as ours, but inserting an imaginary i in the exponent and keeping only a finite number M 
of nonzero arguments 12 We use the shorthand 2^ M) (/i,jt/, T) when there is no possibility of confusion. As we shall 
see,the limit L — > oo, e — > of this object leads to a hierarchy of equations connecting different values of M. We 



sketch the derivation of these equations, with more details in Appendix B and then analyze their basic properties 



We shall not attempt to identify here the precise topologies on the spaces of A's and n's that would permit a rigorous formulation and proof of 
our results. 

1 1 It should be noted that there may be non-realizable solutions of eq. |94[ that do not have this form. In general, statistical moment equations may 
have "parasitic" solutions that do not correspond to realizable solutions of the underlying statistical problem. For an example of this phenomenon 
in the Kraichnan passive scalar model, see 1321 . 

12 Nazarenko employed a different definition in 1131 than that used earlier in (5]H) because he realized the latter definition to be inadequate to 
derive the 1-mode PDF equation (private communication). 
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3.2.1. Derivation 

The method is the same as before. We use the perturbation expansion in e giving the formula ( |75] > for the generating 
functions, with the definitions ( |76| ) of Xl(A,^i) and f77]i-(|8"T) of the JT's. As a consequence of \\ 1 1} all separate terms 
in the prefactors of the JT's are of the same order: 



(112) 



^L(^i-1) + 

K 2 



-^7 — , /ik,/tk 2 , (/Ik, + ~T~, )~j - CA^ )• 



(113) 



Hence, we only need to calculate the leading-order graphical contributions. We already analyzed these in the previous 
section and this discussion carries over here, except that now some wavenumbers are discrete and take on only M 
values (mode 1 for J~\ - J4 and modes 1,2 for SFs) whereas all others are continuous in the infinite-box limit. One 
must consider carefully whether free wavenumbers in graphical su mmations are d iscrete or continuous to see whether 
their contribution is O(M) or 0(L d ). This analysis is carried out in Appendix B It is found that X^(A,fi, T) in the 
large-box limit gets no contributions from J7i, J4 and 0(1) contributions from j 2 > ^3,^/5- The result is 

X m (A,li,T)-X m (A,ii,Q) ~ -9^iS ufi J] J] J d d k 2 d% ^(a ■ kp\h T (a • cofk^f 

j=\ CT=(- 1,0-2,(T3) 

a 



+e 



X 
M 



ft ffL -r(M+2) 



18/5, 



dA kj 3A 2 dA 3 u 3 =^ 3 =o 
>-°Z Z °" 2 f d%d%6 d (o Z -k j )E T (0,gi- co{k.)) 



;'=1 er=(-l,a-2,a- 3 ) 

a2-r(M+l). 

x\Hr\ A k — 

1 1 dA 3 dAj 1 



+ (2^3) 



j=1 m*j,-j £=(-1,0-2,0-3) ^ ^ 



x£ 7 -(w(k j ) + w(-ky), £ • w(k ,))#" ff^^ / a. J„ 



Ai 




^-'i } (2 ~ 3) 



M / a 1 \ r - - 

+ 186 2 J] ^,1(^,1 n ^ ~ 'T" ~ 7TT~ Z d d k 2 d%6 d (g: • k .) 

y=l m*l,-l ' 1 1 ' £=(-0-2,0-3) ^ 

xA r ((r • u(k.))A T (-o- ■ <o(k'))H-<^H + ^fl I JlJ^jJ^" A ' J ") . (1 14) 

* — —J kj.k2.k3 -k,,k 2 ,k3 \ ' J 



Here k^ = (kj,k 2 ,k 3 ) and = (— ky, k2, ks). Taking the small-e limit using the asymptotic formulas (90i, one finds 
the terms proportional to S^o contain secular contributions while the terms proportional to S^iSp . 1 are nonsecular. 
Using ^ to calculate Z (M \T) - Z (M) (0), replacing (Z (M) (T) - Z (M) (0))/T by Z [M) and using time variable r = e 2 f, 
one finally obtains 



d 
dr 



m „ 

Z (M) (A,fi, t) = -36ot^,o 2 I d d k 2 d% 6 d (gz ■ kpS(a ■ <u(kp)|#£ 



j=\ ^=(-1,0-2,0-3) 

d d 2 z {M+2) 



(Ay + ^7-— ) _ _ _ _ -cr 2 A,^- 
1 7 9A ; dA 2 dA ? U2=M=o J m. 



dZ (M+l 



3A 3 dAj u '=o 



dZ (M + l) 
- CT 3 Ay— 

(9A 2 dAj 1 



(115) 



Equation ( 1 15 1 is the main result of this section. Its conditions of validity are similar to those of the spectral hierarchy 
discussed in section U. 1.11 
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3.2.2. Properties 

We now consider the important properties of the multi-mode equations (TT5J. As for the spectral characteristic func- 
tional in the previous section, the factors 6^0 imply that the RP property is preserved in time. Therefore, we can 
consider the generating functions for the amplitudes alone, obtained by setting fi - 0. These form a hierarchy of 
equations, for M = 1, 2, 3, .... 



d M C 2 

— Z (M) (A t) = -36/tt ^ ^ J d%d% 8 d {a ■ kj)S(a ■ 6>(kp)|fljjpJ 



j=i tr=(- 1,0-2,0-3)' (116) 

(/I; + A,- — ) — - — - — _ _ -<T 2 Aj— _ -o- } A;— _ . 

1 ] dAj dA 2 dA^^-=A 3 =o dA 3 dAj^=° dA 2 dA } ^=°> 

By straightforward Fourier transformation in the A variables, one can obtain an equivalenj^jhierarchy of equations for 
the joint PDF's 'P (M, (s 1 , s M ; ki, k M ). For each M - 1, 2, 3, ... these are equations for conservation of probability 

M o 

a 



, ds„ 

m=\ 



ptm + y o r m = 0j (117) 

with the probability flux 

T,i M) = -36ns m fd%d%6 d (g 1 -k m )6(g : -c l j(^ n ))\H^ 2 

r Q>p(M+2) 

ds 2 dsT,— (s u s M , s 2 , s 3 )s 2 S3 

L J os m 

+cr 2 j ds 3 P^ M+1 \su.:,s M ,S3)'s3+o- 3 J dS 2 r (M+l \s x , s M ,7s 2 )s 2 ]. (118) 

We obtain no closed equations for with any choice of M but instead an infinite hierarchy, in which the evolution 
equation for 'P (M) contains and <p( M+2 \ for M = 1,2,3, .... Our results are therefore not in agreement with 

previously proposed closed equations for !P (2) with M - 2 [23 p] 

We do, however, recover the equations for V m and Z, m with M — 1 which were previously obtained 0|6l|7l, 
under appropriate conditions. Namely, assume that the initial data for the hierarchy \\ \6\ are factorized 

M 

Z {M \A U ..., A M , 0; kj , k M ) = \~[ Z{A„„0; k,„), (1 19) 



as would follow from RPA initial conditions, for example. It is then easy to show by substitution into ( 1 16 1 that there 
are solutions which remain factorized 



u 



Z {M) (A l ,...,A M ,T;k l ,...,k M ) = Y]z(A m ,T;k m ), t > 0, (120) 



m=l 

where the factors Z(A, t; k) satisfy the closed equations 



d -Z(A;k) = i^A\l+A—\z(A-k)-y^A—(A-k) (121) 



13 Although naively equivalent, the formulation in terms of amplitude s variables may be preferable to the /i-variable formulation if the ideas 
of (51 I15I are correct that the equations become invalid at sufficiently large amplitudes. In that case, only the .s-formulation would permit the 
incorporation of suitable boundary conditions. See section [4~Tj 

14 Their 2-mode equations (26)-(33) reduce to the asymptotically exact 1-mode equation if one neglects some additional terms that vanish in the 
limit L — > 00 and if one assumes that f^ 1 = yOyW. In this limit in which the 2-mode equation of |23] becomes exact, therefore, it contains no 
more information than the 1 -mode equation. 
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with 



77k = 36n ( d d k 2 d d h <?V ■ k)£(iT • w(k))|i/^ 2 n(k 2 )n(k 3 ) > 0, 

cr=(-l,o-2,o- 3 ) 

y k = 36tt J] f A 2 A 3 <5V ' k)<5(E ' w(k))|^| 2 [tr 3 n(k 2 ) + cr 2 n(k 3 )]. 



(122) 



(123) 



(F— ( — 1 ,0*2,<T 3) ' 



These results are equivalent to the existence of solutions of the PDF hierarchy ( 1 17 1,( 1 18 1 that remain factorized 

M 



r (M) (s u s M , t; kj , k M ) = Y] P(s m , r; k m ), r > 
where the 1-mode pdfs P(s, t; k) satisfy 



m=\ 



8t 



d_ 
ds 



dP 



(124) 



(125) 



Preservation of RPA therefore follows if one can prove uniqueness of solutions of the hierarchy \\\6\ or, equivalently, 
( 1 17 1, ( 1 18 1. The equations for the factors agree with previous results for the 1-mode equations ||5]|6|I3- 

The equations ( 125 1 are not simple linear Fokker-Planck equations, however, but are instead nonlinear Markov 
evolution equations in the sense of McKean [25]. That is, the solutions must satisfy a set of self-consistency condi- 
tions, 

n(k,r) = j ds sP(s, t; k) (126) 

where n(k, r) is the same spectrum that appears in the formulas for the coefficients (122 1,(123 I. These equations 
are the exact solutions of a model of "self-consistent Langevin equations", like those for the DIA turbulence closure 
11331 l34l . Here the model equations take the form of the stochastic differential equations 



ds k = (77k - 7kSk)dT + ^2rj k s k dW k , 



(127) 



interpreted in the Ito sense and with self-consistent determination of n(k, r) via (126i. Solutions of this stochastic 



(128) 



model can be realized by a Monte Carlo procedure of McKean 8251, with ( 127 1 generalized to 

ds k n) = (7?k - y k 4"Vf + yj^sf dW k "\ k€A>=l, —,N, 
with the spectrum obtained by an A^-sample average 



1 N 

1 V («) 



(129) 



and with y k , rj k given by formulas ( |122| i,( [T2"3"j ) in which wavenumber integrals are discretized as sums. Taking first 
>> 1 and then L » a yields a solution of the PDF equations ( 125 1. This procedure works as well to solve the 
general hierarchy equations ( 1 17 1,( 1 18 1, without assuming factorized initial data. For the factorized case, a far simpler 
procedure is to solve first the wave kinetic equation Q for n(k, r) and then, using this as input, to solve the 1-mode 
equation ( 125 1 in order to obtain P{s, t; k) for any wavenumber mode k of interest. As initial condition one may take 
any P(s, 0; k) which satisfies the consistency condition ( |126| l at r = 0, since this condition is preserved in time by the 
1 -mode equation. 

It is interesting that, independent of the initial condition P(s, 0; k), the solution P(s, t; k) relaxes as t increases to 
a Rayleigh distribution 

G(*,r,k)= — exp(-s/B(k,T)), (130) 
n(k, t) 

which corresponds to a Gaussian distribution of the Fourier coefficient 7i k (T). It is easy to check that, for any solution 
n(k, t) of the wave kinetic equation Q, Q(s, r;k) defined above solves the 1-mode PDF equation (125 1. (Since 
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wavenumber k appears only as a parameter in our argument, we suppress its appearance for the rest of this paragraph.) 
The relaxation of a general solution P to Q is indicated by an //-theorem for the relative entropyY^\ 

H(P\Q) = j ds P(i)ln^^j = j ds P(s) In P{s) + Inn + 1. (131) 

This is a convex function of P, non-negative, and vanishing only for P — Q |35|. Taking the time-derivative using 
( |125| i, it is straightforward to derive 

d vrvr wnr w s\8 s P(s,t)\ 2 rj 

—H(P(t)\Q(t)) = -r) ds — — — + — -, (132) 

dr J P(s,t) n(T) 

where 

J -sd s P(s,r)ds = jp(s,T)ds=l (133) 

was used to cancel terms involving the y coefficient. But note the self-consistency condition n(r) = J sP(s,T)ds 
implies 



^-H(P(t)\Q(t)) = -tj 



s\d s P(s, t)| 2 1 



P(S,T) j sP(S,T)ds^ 

The inequality follows from the Cauchy-Schwartz inequality applied to ( | 1 33 1 



< 0. (134) 



/ ^ ' ^ ( ~ dsP) ds - "vj sP(s) ds ' J S ^P~ ds - (135) 



Equality holds and relative entropy production vanishes if and only if yfsP = c Jj(—d s P) for some constant c, or 
P = -cd s P. The solution of this latter equation gives P — Q with n = 1 /c. We thus see that P(r) should relax to Q(t) 
as t increases. For the purpose of later discussion, we emphasize that this argument assumes that the kinetic theory is 
valid over the entire range of amplitudes s e (0, oo) and it could otherwise fail. 

Let us remark that the relative entropy has a simple probabilistic meaning, similar to that discovered by Boltzmann 



for his entropy function, which involves the empirical PDF Pl(s; A) defined in ( 59 ). For a single cell A in wavenumber 
space A*, let fij^n^ be the microcanonical measure on the set of microstates r^ nr] = {(J,i[t) : [«l(A) - w(A)| < rj}. 
Then, with respect to this microcanonical measure, the most probable value of Pl(s; A) as first L — > 0, then rj — > is 
Q(s; A) = e _s /"(A) /re(A). However, the probability of observing another PDF P as a rare fluctuation is 

lim lim -J— \n^, hn ({P L * P}) = -H(P\Q). (136) 

,_»0 L^ca N L (A) ' 

This can be shown heuristically by adapting the original argument of Boltzmann [30). Let P(s) be a probability 
density function over possible values of s, let S = {Lj, j = 1, ...,R} be a finite partition of the positive reals into 

intervals with Sj the midpoint values and define the integers Nj = N JT P(s)ds where [ J denotes integer part and 

= Nl(A). (Since A is fixed in this argument, we hereafter omit explicit reference to that quantity in the remainder 
of this paragraph.) Note that Nj — N for large enough L. Now let Nj denote the number of modes k e A^ n A 
such that 7k £ £j and define the set Tpx,l = {(J, <p) ■ Nj = Nj, j = 1, ...,R}. The Liouville measure of this set is 
\Tp,i.,l\ = N f': N]jl \I.i\ Nl ■ ■ ■ \Lr\ Nr ~ exp{-NZj \Zj\PjlnPj} with P } = Nj/N]£j\ and using Stirling's approximation for 
the factorials, exactly as in the original argument of Boltzmann. Then 

lim lim -lnirpvd = - f ds P(s)iaP(s). (137) 

|£[->0 h» ' ' J 



15 A related set of observations were made in | 7 |, section 9. 
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Now consider the set F„ 



nji,L - [ ~ > {p-\n-Y ll \Zj\sjPj\<ri^'pX,L- I ts Liouville measure |r„. ); /J is dominated in the limits L — > oo, 
r] — > 0, |S| — > by |rp£ Ll for the pdf P satisfying J ds sP(s) - n with the largest entropy - J ds P(s) In P(s). As is 
well-known, this maximum entropy distribution is the Rayleigh pdf Q(s) - e~ s ^ n /n. Thus, 



lim lim lim 

|Z|-»0 L-*ca 



in|r„, 



I 



ds QO) In gO) = Inn + 1. 



(138) 



■(/ dsP(s)lnP(s) + ln« + l) , as claimed. 
The result (1136b becomes the standard Sanov theorem [35 1, if the measure fi&,„,r, is replaced by the RPA measure for 



Taking n„J{P L * P}) 



\r„.„.t\ ■ 



then limpi^o lim^o liniz,- 



n 111 [r„,i| 



which each mode J^, k e n A has the independent density Q(s; A). 

A simple picture emerges of entropy production in kinetic wave turbulence. For the macroscopic state specified 
by the 1-mode PDF's P(s; k) of all wavenumbers k, the Boltzmann entropy obtained by counting microstates is 



S(P) 



jds j d d 



k P(s\K)\nP(s\\i), 



which satisfies the 2nd-law dS /dr > under the evolution equations ( 125 i. Associated to the decomposition 



S(P) = -H(P\Q) + S(Q), 



(139) 



(140) 



there are two distinct sources of entropy production. The increase in the usual entropy of wave turbulence S{Q) = 
J d d k (ln«(k) +1) corresponds to the approach of the wave spectrum n(k) to its thermal equilibrium form n eq (k) = 
r/w(k). The increase in -H(P\Q) corresponds to the approach of the PDF's P(s;k) to the Rayleigh distributions 
Q(s; k) with the same spectrum n(k) - J ds sP(s; k). It is interesting that this approach to Gaussian statistics is 
quite slow according to the equation ( 125| l, taking place on the nonlinear time-scale 0(e 2 ). This is in contrast to the 
physical space statistics which — even if non-Gaussian initially — are expected to approach Gaussian on the fast wave 
time-scale. See footnote ^[T] The physical-space statistics are furthermore expected to remain close to Gaussian under 
the influence of the weakly nonlinear dynamics. There is no contradiction with our //-theorem, since preservation of 
RPA in the weak senses (57 1, (62 1 and the results in section [2] imply that the physical space statistics of the random 
field m(x, t) indeed remain Gaussian to leading orders in e and 1 /L Il24l . 

The above argument motivates a microcanonical measure //p,a,l,/;,l defined for a given P(s;k) by normalizing 
Liouville measure on the set 



^p,ax,i,l = i (J, <A) 



I dsP L (s, A,-)- ( dsP(s,Ai) 
Js, Jz, 



< ij, i= l,...,P,j = l,...,R 



(141) 



for partitions A of A* and £ of M + . This measure satisfies the RP property, as it contains no dependence on the phases. 
One can define for this measure the M-mode correlation functions of the empirical PDF, Pt(<^) = EkeA^ V(A, k): 

.., <Pm) = fiPA£, n ,L (Pl(<Pi) ■ ■ ■ Pl{^m)) (142) 
(Pm on M. + x A* vanishing at infinity. The arguments of ETI can be 



r (M> (wi 

r >,A,£,/7,£ W1 ' 



for any choice of continuous functions ip\ 
adapted to show that 



lim lim lim V^} t„ l (<Pu -,<Pm) = P(<Pi) • • • P(<Pm) 



|A|,|£| 



(143) 



where P(<p) — C ds f d d k tp(s, k)P(s; k). Thus the factorization property of the M-mode PDFs holds asymptotically 
for the microcanonical measure fip^x^.L, although it does not satisfy the strict RPA property. Our derivation of the 



PDF hierarchy equations ( 1 17 1,( 1 18 1 assumed existence of the M-mode PDF's in a somewhat stronger sense in the 
limit as L — > oo (i.e. for fixed values of ki, ..,k M ), but it is reasonable to expect that they remain valid for such 
"generalized RPA" initial conditions. 

Lastly, we can ask what are the solutions of the PDF hierarchy for general initial data which are RP but which are 
not factorized even asymptotically. One must assume at least that 
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dsi--- ds M d%---d d k M 
o Jo J 



X ^iCsi.ki) • • ■ <Pm(sm, k M )P { M) (s U -, «m; k! , .., k M ), (144) 



in order to provide suitable initial data P^{s\ , .., Sm\ ki, .., k^) for the hierarchy. However, with appropriate technical 
assumptions, it follows by the argument of Spohn [22 1 that 

^ M) (ii,..,5M;ki ) ..,k M ) = fdp (P)P(s 1 ;k 1 )---P(s M ;k M ), (145) 



where po is a probability measure on the PDF's. Since the PDF hierarchy equations (117i,(118i are linear in the 
!P (M) 's, a solution will be provided by 



^ M) (a 1 ,..,a J ,,T;ki,..,k Jf ) = J d Po (P) P(s U T;kr)---P(s M ,T;k M ). 



(146) 



This will be the only solution if existence and uniqueness of solutions holds for the PDF hierarchy. Thus, the most 
general realizable solutions of the PDF hierarchy equations are again expected to be "super-statistical solutions" that 



correspond to ensembles of solutions of the 1-mode PDF equations ( 125 1 with random initial conditions 



4. Intermittency in Kinetic Wave Turbulence 

One of the important potential applications of multimode equations in wave kinetics is the explanation of observed 
intermittency and anomalous scaling in wave turbulence lfT6l[T5l[T7ll . On the face of it, wave kinetic theory appears 
to have few resources to explain such phenomena and instead has all the attributes of a "mean-field theory" ll36l . As 
we have seen in previous sections, it is a theory which ignores fluctuations and in which important quantities in fact 
become deterministic. All interactions of the infinite collection of modes is through a mean-field n(k, r). Furthermore, 
scaling exponents in the Kolmogorov solutions of the wave kinetic equations are given by dimensional analysis [ 18 1. 
These are the hallmarks of a mean-field theory, which generally cannot predict anomalous scaling. 

Nevertheless, there are at least two approaches based on wave kinetics which seem to have some promise to explain 
intermittency and non-Gaussian statistics. One is the idea of a "cascade in amplitude space" that was proposed by Choi 
et al. ll5l [T5ll . based on the 1-mode PDF equation ( 125 1. Another is the idea of a "super-turbulence" in wave-kinetics, 



advanced in the present work. We shall now discuss both of these possibilities. 
4.1. Cascade in Amplitude Space? 



As already mentioned at the end of section 3.1.1 wave kinetic equations generally cannot be uniformly valid over 
the whole range of wavenumbers, but must break down in either low or high wavenumbers where nonlinearity become 
strong [ 18 1. Choi et al. (also llT5l ) have proposed that there is a similar non-uniformity in amplitude space, with 



equation (125 1 for the 1-mode PDF P(s,t;V) restricted in validity to s < s„i, where s„i is the amplitude for which 
nonlinear interactions become strong at wavenumber k. [For the remainder of this section we shall assume that the 
wavenumber k is fixed and drop it as an explicit label.] For amplitudes s > s„i strong-interaction processes such 
as "wave-breaking", "cresting", etc. are assumed to occur which are beyond the description of wave kinetics. The 
hope is that the effects of such nonlinear processes can be modeled by supplementing the equation ( 125| l with suitable 



boundary conditions and forcing terms. We find this a very intriguing suggestion but, as we now argue, it seems to us 
to have as yet no successful formulation, either analytically or physically. 



Let us review the specific proposals of 1151 [151 . They note that eq.( 125 1 is a probability conservation law 

d,P + d s F = (147) 

with probability flux in amplitude space given by 



F(s) = -s(yP + rjd s P). 
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(148) 



The general time-independent solution of ( 147 i with constant flux F (s) = F* was observed in [5 1 to be 



P(s) = Ce- S/V - —Ei(s/v) e - s/v (149) 

with v = rj/y and with Ei(x) = P f* e' j the standard exponential integral. The first term has zero flux F(s) = while 
the second term has non-vanishing flux F(s) = F„. Since Ei(x) > for x > x* = 0.3725, positivity of P(s) for s » v 
requires F* < 0. When F t = then = Q(s), the Rayleigh distribution with spectral density n — v. Not noted in 
is the fact that one must more generally make a distinction between v and n given by the self-consistency condition 



ds sP(s), (150) 



which are only the same when P = Q. The above integral diverges for F* + 0, if the upper limit extends to infinity, 
because P(s) ~ due to the asymptotics Ei(x) ~ ^(l + ^r + f?+""") f° r x » 1. The integral is finite if it is cut off 
at the upper limit s„i. This assumes that whatever nonlinear processes such as "wave-breaking" occur at s > s„i lead 
to a PDF which decays very rapidly for s > s„i. 

The papers If5l [l5ll propose to model the nonlinear effects at s > s„i by a negative flux F(s„i) = F , < 0, which 
represents a flow of samples back into the kinetic range s < s„i, at least in some range of wavenumbers. For example, 
this might occur at low-wavenumbers as wavecaps and cusps in the nonlinear range break up and feed back into 



the weak, incoherent background. However, the equation ( 147 1 will then no longer conserve probability but instead 
satisfies j t J Q "' dsP(s) = \F*\ > 0. Paper [5] proposes to add to the equation a constant "drag" D 

d t P + d s F = -D, s<s n i (151) 

to represent the dilution in the weight of samples at s < s„i by the flux F* of new samples into the ensemble at s = s n ). 
It is easy to see that one must choose D = \F*\/s„i to conserve total probability^ Paper then claims that the general 
solution to the modified equation is P(s) = Ce~ s ^ n + ^Ei(s/n - log(s/nJ)e~ s ^ n . The picture proposed in |j5][T5] is 
one of an "inverse probability cascade" with flux F* < in amplitude space, coexisting with the usual forward 
energy cascade in wavenumber space. The power-law tails P(s) ~ ^ for s » v are a source of non-Gaussianity and 
intermittency in wave turbulence if F, + 0. A "critical balance" argument is used in 1151 to estimate that |F»| ~ yn/s„i. 

It is easy to see, however, that the proposed cascade picture cannot be correct for the specific model in 0, on both 
physical and mathematical grounds. A constant "drag" D is not localized at small and large s but has effects felt over 
the whole range of s. Thus, one would not expect a constant-flux "inertial range" to exist for such a drag. In fact, for 
any choice of D, a particular solution of the inhomogeneous equation ( 151 1 is 



P parti*) = D/y. (152) 



The general solution P(s) of ( 151 1 is a superposition of P par t(s) with any solution of the homogeneous equation ( 147 1 



P(s) = ° + ^± e -'/r + %-s/v m { S \, (153) 
J V V \V I 

This has flux F{s) = -Ds - Ciy. With the b.c. F(s„i) = F* and the unique probability-conserving choice D = \F m \/s„i 
one finds Ci — and thus 

using normalization. Note the requirement |F»] < y for positivity of this solution, which is then a superposition of a 
Rayleigh distribution and a uniform distribution. As expected, it has not constant flux but instead F(s) = F t s/s„[. 
We claim that (154i is the correct solution of the model formulated in Q. It gives an alternative possibility to 



explain intermittency and non-Rayleigh distributions, with constant tails P(s) oc s for s » v rather than P(s) oc s 



16 In fact, 1 5 1 proposed to take D = -yP(s n f)/ s„i, which does not conserve probability and which is not even dimensionally correct. This equation 
a typographical error not caught by the authors of | 5 1 (S. Nazarenko, private communication). 
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Using the estimate |F»| ~ yn/s„i of [ 15] one finds that the realizability inequality < y is satisfied whenever n < s n i. 
However, it is unclear to us that the flux b.c. and constant drag D is a physically correct model of wave-breaking. 
Constant-flux solutions as originally proposed in [5] could be relevant with a more general model of strong nonlinear 
effects, in which the constant D is replaced with a function D(s) such that "' ds D(s) = \F t \ but for which D(s) nearly 
vanishes except at the upper and lower limits. This would provide a "transparency window" that would allow constant 
probability-flux solutions to exist. 

Whatever may be the correct model of the strong nonlinear effects at s > s„i, an important general point is that 
any modification of the 1-mode PDF e q.fll4 7|> requires that the wave kinetic equation must also be modified. This can 
be illustrated for the specific model in ( 151 1 with D = \F,\/s„i. Taking the time-derivative of n = J ds sP(s) gives 



-yn + i] - Tjs nl P(s n i) + -\FJ\Sni, 



(155) 



which differs from the standard kinetic equation h = —y n + rj. The stationary form of the "modified kinetic equation" 
( 155 1 is easily checked to be valid for our solution ( 154 1 by calculating n — L "' ds sP(s) to be 



V-VS„lP(S n i) + —S n l 

2y 

V| 1 + — S n l, 

7 ) 2y 



(156) 



with the latter approximation valid for s„i » v. The correction to the usual value n — v predicted by the kinetic 
equation need not be small. If we use the estimate ~ yn/ s m i from "critical balance" [ 1 5 1 then the second term in 



( 156 1 is comparable to the first term, or even much larger if v «: n. If such changes in «(k) occur for a large set of 



wavenumbers k, then the coefficients y and rj defined through the integrals ( 122 1,( 123 1 may also be strongly modified. 
These remarks make clear the nontriviality of constructing a self-consistent hybrid model of wave kinetics and of the 
strong nonlinear effects. 



4.2. Super-Turbulence of Wave Kinetics? 

A second possibility to explain intermittency and anomalous scaling entirely within the framework of wave kinet- 
ics is by "super- turbulence". We shall discuss this mechanism here, mainly at a general, qualitative level. 

It is useful to begin with a review of the work of Spohn [22], who already discussed turbulence in gas dynamics as 
a natural motivation for ensembles of Boltzmann solutions. This is easiest to understand in the hydrodynamic regime 
where length-scale of variations in local equilibrium parameters (temperature, density, velocity) is much larger than 
the mean-free-path length. The Boltzmann equation reduces then to a hydrodynamic description and one can appeal 
to the extensive literature on turbulent solutions of hydrodynamic equations. The specific example discussed in ll22l is 
the Rayleigh-Bernard system, considered as a many-particle system subject to thermal boundary conditions. In such 
a situation driven by the boundary conditions, Boltzmann's original //-theorem is no longer valid and the stationary 
solution of the Boltzmann equation for the pure conducting state is no longer purely Maxwellian but has a small 
correction corresponding to thermal non-equilibrium [37|. In the turbulent regime at sufficiently high Reynolds and 
Rayleigh numbers, the laminar purely conducting state is unstable and turbulent convection develops. The temporal 
dynamics is chaotic so that long-time-averages, for example, are described naturally by invariant measures corre- 
spondng to ensembles of Boltzmann solutions. 

Such a statistical description is natural also for turbulent situations without driving by body forces or boundaries 
and with time-dependent statistics. For example, a single realization of decaying, statistically homogeneous turbu- 
lence is very spatially complex and heterogeneous. Averages over large volumes — by space-ergodicity — are again 
described by ensembles of solutions. Note in this example that Boltzmann's //-theorem applies, with a monotonic 
increase of entropy due to heating of the fluid. There is eventual approach to global thermodynamic equilibrium and 
a space-independent Maxwellian 1 -particle distribution, but the turbulent state with strong spatial variations exists 
as a long transient for an intermediate range of times. The previous examples are in the hydrodynamic regime, but 
there should be similar turbulence in the fully kinetic regime. For example, in compressible turbulence in gases at 
sufficiently high Mach numbers the thickness of shocks should be of the order of the mean-free-path length |38| and, 
for a quantitative description, the Boltzmann equation should be used rather than hydrodynamic equations. 
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These ideas can carry over from the Boltzmann equation to the wave kinetic equation. As noted in the introduction, 
Zakharov et al. ([ Q, section 4.2.2) had already pointed out that strong instability of the Kolmogorov solutions to the 
wave kinetic equation could lead to such a "secondary turbulence" . The statistical distribution of wave amplitudes that 
would be obtained in a long time-series from experiment or simulation would then correspond to a random ensemble 



of solutions of the 1-mode PDF equation ( 125 i: 



!P (1, (s;k) = j dp(P) P(s;k). (157) 

Even if the individual solutions P(s; k) were close to Rayleigh distributions, the measurable distributions !P (1) (s; k) 
could be arbitrarily far from Rayleigh, depending upon the "super-statistical" measure p. To determine this measure 
becomes a very difficult problem, as daunting as the corresponding problem for Navier-Stokes turbulence! Here we 
can only identify some of the sources of randomness, including instability, random forcing, and boundary conditions. 
We consider these in turn. 

Instability of the Kolmogorov solution is a natural origin for "super-turbulence". As a related example, consider 
the GOY shell model of turbulence, a dynamical system with a stationary Kolmogorov solution supporting a constant 
energy flux with dimensional scaling [39]. This solution, however, is subject to a rich array of instabilities |j40ll4"Tl 
and the statistical behavior of the GOY model shows strong intermittency and anomalous scaling. There is a detailed 
theory of linear stability of the Kolmogorov solutions of wave kinetic theory P2l |43l l44l l45l . reviewed in [1 1, Ch. 
4. This theory studies perturbations to a scale-homogeneous Kolmogorov solution which, strictly speaking, requires 
an infinitely long inertial interval. The conclusion of this theory is that, for most common cases, the Kolmogorov 
solution is linearly stable. However, the theory, although sophisticated, is not quite definitive. Quoting from Q: 

"It should be noted that a more consistent formulation of the problem of the stability of the Kolmogorov 
spectrum should be as follows. First of all, the kinetic equation should be supplemented by terms de- 
scribing the isotropic pumping and damping regions and a stationary solution of this equation should be 
found that is close to the Kolmogorov spectrum in an interval (fci,^); outside this interval the solution 
may strongly differ from the Kolmogorov spectrum. Then the kinetic equation must be linearized in the 
vicinity of this stationary solution and expanded in angular harmonics. ■ • ■ Having examined the behav- 
ior of the solutions equation, one should clarify the changes that occur when the ranges of the source 
and sink in A:-space go to zero or to infinity and examine the behavior of the perturbations established in 
the interval C&i.&i). Finally, one should analyze in which situations this behavior is independent of the 
specific type of the source and sink. 

This program for examining the stability of Kolmogorov spectra turns out to be too complex. Currently 
there exists no strict proof of the fact that in general the kinetic equation with a source and sink has a 
stationary solution close to the Kolmogorov spectrum in some interval." 

In fact, the situation is even less certain than this statement implies. It has recently been found that the space- 
independent Kolmogorov-Zakharov solutions can be linearly unstable to spatially inhomogeneous perturbations [46|. 
In addition, there are hydrodynamic flows such as plane Couette and pipe Poiseuille which, despite being linearly 
stable for all Reynolds numbers, become turbulent in laboratory experiments and in numerical simulations at moderate 
Reynolds number 1371 . The usual understanding is that these flows are unstable to finite amplitude perturbations, e.g. 
|48|. Thus, linear stability does not rule out transition to turbulence. For these reasons, we must regard "super- 
turbulence" as a viable possibility in wave kinetics. 

As emphasized in the previous quote, external forcings are required to maintain a stationary energy cascade state. 
It is important to consider more deeply the origin and role of such forces. The monograph Q considers the necessary 
conditions on the force for the existence of a stationary distribution (section 2.2.3) and how to match the Kolmogorov 
solutions to the wavenumber regions of the force (section 3.4). This discussion assumes a particular idealized model 
for the forcing, in which to the collision integral /(k, t) of the wave kinetic equation there is added an additional term: 

n(k, t) = F(k, T)n(k, r) + /(k, r). (158) 

This model describes "pumping" for wavenumbers at which F(k, t) > and "damping" for F(k, f) < 0. Validity of 
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this model requires a suitably weak, slowly-changing force in the equations of motion (31 1: 



(159) 



Here • • ■ denotes the original terms in (31 1. The addition of this term to the dynamics can be easily accommodated 



into our derivations, with a new term appearing in eq.(73 1 



in eq.(79i 



in eq.(94i 



dZ 

(It 



J3 = 1%^,^;, + •• 



' J k T(k, T)A(k) 



6Z 
6A{k) 



in eq.(115l 



dT 



8A„ 



(160) 



(161) 



(162) 



(163) 



and in eq.( | 1 1 8| > 



ri M) = r(k m ,rK + 



(164) 

Thus, the standard forcing model from [ 1 ] is obtained. All of our previous results, e.g. on preservation of factorized 
solutions, law of large numbers, etc., still carry over. In general, forced wave turbulence can be described by kinetic 
equations only for suitable assumptions on the forcing, like those above. In the above derivation, the only randomness 
was in the phases and amplitudes of the waves. However, the forcing amplitudes F(k, t) may be chosen to be a 
realization of a stationary and time-ergodic random process. In that case, long-time averages in the statistical steady- 
state can be described by averages over ensembles of solutions of the wave kinetic equations with different realizations 
of the force. In fact, [49 1 find that in their experiments the forcing is well modelled by an Ornstein-Uhlenbeck process 
with large fluctuations in the energy input rate. Even without chaos or turbulence in the wave kinetic equations, this 
will yield "super-randomness", since the Kolmogorov solution is neutrally stable to change in the energy flux rate (TJ. 

In addition to pumping by an external force, hydrodynamic turbulence may also be (in fact, more commonly, is) 
driven by boundary conditions. Such a situation can also occur in wave kinetics, as we have discussed in the previous 
section 4. 1 where strongly nonlinear effects such as wave-breaking can provide both sinks and sources to kinetic 
wave turbulence. These effects restrict the validity of wave kinetics to a "bounded domain" in wavenumbers and 
amplitudes and could provide suitable forcings at the boundaries to represent those nonlinear effects. In this setting, 
Kolmogorov cascade solutions which are stable when considered without restriction, might become unstable and 
"super-turbulence" develop. The two scenarios in this section and the previous one are thus not completely exclusive. 



5. Conclusions 

The major contributions of this paper are the following: 

1 . We have derived by formal calculations the leading-order, asymptotically correct multimode equations for wave 
kinetics in Hamiltonian systems with 3-wave resonances. 

2. We have shown that these equations possess factorized solutions for factorized initial conditions, corresponding 
to "propagation of chaos" or preservation of "random phases & amplitudes". 

3. We introduced the "empirical spectrum" and "empirical 1-mode PDF" and showed that the above factorization 
implies a law of large numbers, so that these quantities satisfy the wave-kinetic closure equations for nearly 
every initial realization of random phases and amplitudes, without necessity of averaging. 

4. We have demonstrated the close formal relations of wave kinetic theory with the kinetic theory of gases, espe- 
cially regarding the role of entropy and the 2nd law of thermodynamics in both theories. 
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5. We have completely classified all realizable solutions of our multimode equations and have shown that they cor- 
respond to "super-statistical solutions", or ensembles of the wave-kinetic closure equations with either random 
initial conditions or random forces. 

6. We have exploited our results to discuss the possibilities to explain intermittency and non-Gaussian statistics of 
wave turbulence within the kinetic description, in particular by a "super-turbulence" corresponding to chaotic 
or stochastic solutions of the wave-kinetic closure equations. 

Our work solves the problem posed in [2], section 5.2.5 on the continuum limit of finite-dimensional wave turbulence 
to show "how the natural closure arises in taking the limit L — ♦ oo." Although we have focused here on 3-wave 
resonant models, we expect that analogous results should hold also for other closely related wave systems, e.g. those 
whose first nontrivial resonances are 4-wave. 

It is perhaps important to emphasize what we have not accomplished in this paper: we have not provided a rigorous 
mathematical proof of the wave kinetic theory. At this time, no set of sufficient conditions are known that would imply 
the validity of wave-kinetics for any general class of systems (12], section 5.2.6). Our derivations fail to constitute 
a proof in particular because we have made no attempt to rigorously bound the 0(e y ) remainder term and higher 



order terms in the perturbation expansion (76 1. In principle, there exist methods to calculate the collision integral 
in the wave kinetic equation formally to any order in the small parameter e 11271 l50l I5D . As a matter of fact, there 
are reasons to expect that the collision integral is non-analytic in e, if the kinetic theory of gases is any guide. The 
multiscale asymptotic formalism devised by Bogolyubov [52| to calculate systematic corrections to the Boltzmann 
equation in powers of gas density is closely analogous to the methods employed to derive the wave kinetic equation 
perturb atively in e. However, it was discovered in the mid-1960's that the 4th- and higher-order terms in Bogolyubov's 
density expansion of the collision operator for gases are divergent and the true dependence appears to be non-analytic 
in density [53]. Even for gas kinetics, many fundamental issues regarding higher-order corrections remain unresolved. 
The rigorous works of Lanford GUI I2T1 imply that these higher-order terms are asymptotically negligible in the low- 
density limit for gases, but similar mathematical study of wave kinetics is only just begun. For a recent rigorous proof 
of validity of wave-kinetic predictions in a 4-wave model problem at thermal equilibrium, see Ifl4"l . 

Our goal here was the more directly physical one of formulating the limiting equations of wave kinetics and un- 
derstanding their content. We have found that the equation <|3j originally obtained by Peierls [ 8 ] and many subsequent 
authors |9l [10] [3] [TJ is not the leading-order, asymptotically valid equation in the wave kinetic limit. The Peierls 
equation contains additional terms which vanish as L — > oo and which do not appear in our limit equations. It should 
be stressed that the derivation of the Peierls equation in J5j [SI El an d a l so earlier derivations of Peierls himself |8|, 
Brout & Prigogine [9], Zaslavskii & Sagdeev [10], etc. are not systematic in the limit L — > oo. At present, the Peierls 
equation has no established validity (or even a precise meaning) in the large-volume limit. Although the Peierls equa- 
tion does not have any obvious asymptotic validity, it does serve as a "generating equation'^] from which all of our 
simpler equations can be obtained in the limit L — > oo. It is an interesting issue whether the additional terms retained 
in the Peierls equation describe correctly sub-leading phenomena in powers of l/L. This question shall be addressed 



in our subsequent paper [24|, where we shall derive the limiting equation for the generating functional ( 1 10 1 with all 



N modes retained and use our techniques to evaluate the next-to-leading order corrections in l/L to the asymptotic 



limit. In any case, the virtue of our new equations (92 1 and ( 1 15 i is that they contain all and only the leading-order 
terms and thus allow us to clarify the structure of that limit. We have completely characterized the solutions of the 
limiting hierarchies and shown them to consist of random ensembles of solutions of the wave kinetic equation Q 



and 1-mode PDF equation ( 125) . This should permit a better evaluation of the theory by simulation and experiment. 



In particular, our most subversive conclusion is that the "mean-field" scaling of the Kolmogorov-Zakharov solutions 
may not be the true prediction of kinetic theory for scaling exponents of spectra and higher-order statistics in wave 
turbulence and that wave kinetics may instead allow for intermittency and anomalous scaling. 
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This terminology was suggested to us by S. Nazarenko (private communication). 
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Appendix A. Derivation of the Spectral Hierarchy 



As discussed in the text, we have verified the results of |6] for the J'-terms in eqs.(77 1-(78 I, up to minor correc- 
tions. The contributions to the spectral generating functional are obtained by making the substitutions = {2n/L) d Ji 
and A\ = iAQty) and taking the limit L — > oo. The asymptotics of the various terms can be summarized in the follow- 
ing general rule: each action variable J changes the order by L~ d and each free sum over wavenumbers by L d . We 
illustrate this rule in our detailed calculations below. 



Calculation ofSTi-' Substituting the definition (|70|l of a. , one obtains 



1.2.3 



,o"2k2+a"3k3 



(A.l) 



which can be represented by the following diagram before averaging over phases: 



1 « 




After phase averaging, there are two contributions represented by the following diagrams: 





or explicitly written as 

J\ = ^ (^1 + ^J~) L k'uk^ V^2^ A r(Q'2^ l 2 + 0~3^3 ~ ^l)^k|, t r 2 k 2 +o-3k3^ 1 . 1^2,-0-2^3,-0-3 ]"[ <Wo 



1,2,3 



m* 1,2,3 



1,2 1/1 ' ' m*l,2 



(A.2) 



Note that a prime / on a sum indicates that all wavenumbers must be distinct. Here all of the summations are "pinned." 
Now making the substitutions = (2n/L) d Ji and Ai = iA(ki) into the above and taking the limit L — > 00 gives 



^ = \(t) z z^sv¥^^- w ®^*^ n <w 

^ ' oj=(-l,o-2,o-3) k " 1/1 m*l,2,3 

4(Tf 2 z z'^-^^-^^^n^ 

' ' <r=(-l,2o- 2 ) k' V-'l m*l,2 



(A.3) 



where y! = ip.\,H2) and k' = (k;,k2). Hence J'x is 0(L d ^ 2 ) in the large box limit. In fact, the result easily follows 
from the aforementioned asymptotic rules. The amplitude factor V-A-^2^3 changes the order of Ji by Zr 3< ^ 2 . Since 
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there is no free sum here, the contribution of the term proportional to A\ is at most 0{L M / 2 ). The term proportional 
to hi introduces an additional action variable l/Ji and t hus g ives the leading contribution 0(L~ d ^ 2 ). 



Calculation of Si'- Again substituting the definition (72 i of a^\ one obtains 



4ji 

A r (cr 2 <jJ2 + cr 3 w 3 - u)\)A T {cr^U)4 + cr 5 aj 5 + Wi^k^kz+cr^-k^kt+t^ks 
which can be represented by the following diagram before phase-averaging: 

4v A 



1 





(A.4) 



After phase averaging, Si has three types of contributions. The leading contributions are two type I diagrams with no 
external couplings. We show one such diagram here (the other one is obtained by the exchange 2 <-> 3). 




These contain three distinct wavenumbers and two vertices, but the wavenumber delta-functions at the two vertices 
turn out to give the same constraint. Thus the type I diagrams contain two free wavenumbers corresponding to un- 
constrained sums. Subleading contributions are given by type II diagrams with two wavenumbers externally coupled. 
We show one diagram here with the wavenumbers 2 and 4 coupled to distinct external blobs. The three others can be 
obtained by exchanging 4 <-> 5 and 2 <-> 3. 




These subleading contributions contain only one free wavenumber. The two delta functions 5ki,o- 2 k2+o-3k 3 an d ^-ki,o- 4 k4+o-5k 5 
along with the internal connection 5k,,k 5 5o-3,-o- 5 require <r 2 k2 = -0-4^.4. Therefore, k 2 = -k4 (since kj and 1*4 must be 
distinct) and o~2 = cr 4 in this diagram. In principle, another type II contribution could exist with two wavenumbers 
coupled to the same external blob, as follows: 
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together with diagrams resulting from 4 <-> 5 and 2 «-» 3. However, the condition cr 2 k2 = 
that k 2 = k4, cr 2 = — o"4, and fi 2 = 0, which coincides with the type I diagram. 

All other contributions are type III where all wavenumbers are pinned, for example, 



-cr 4 k4 in this case requires 




Now consider large L. The amplitude factor ^|7 2 JyJ~^Jl changes the order by L 2d . Now consider the contributions 
proportional to the three terms in the prefactor A\ + A 2 Ji — In the term proportional to A\ the type I diagrams 
provide two free sums that increase the order by L 2d , giving an 0(1) contribution. The next order contributions from 
type II diagrams have only one free wavenumber and thus are at most 0(L d ). The term proportional to A^Jj is at 

most 0(L~ d ) because of the additional action variable J\ . For the term proportional to jj-, we note for both type I and 
type II diagrams k\ is not pinned, which requires \i\ to be zero. Thus the only nonzero contribution here must come 
from type III graphs, which is at most 0(L~ d ). Therefore we can conclude that the leading contributions in £Ti are 
0(1), or explicitly 

Ji ~ 96^Y j A 1 J 2 J3\^ 2 2\ 1 \At(o-20J2 + o-3CJ3-oji)\ 2 (A.5) 

1,2,3 



Substituting = Cl-KjLfj; and A\ = z'/l(ki) and taking the limit L — > oo using (84i, one finds 

I d"k l d"k 2 d"k 3 d"(a ■ k)|A r (cr • w(k)j' 2 

0"=( — 1,^2^) * 



( e 2Wj- 2 ) 7 ~ -9/^,o J] f d%d d k 2 d%6?(a-^h T {a-a)Q$)f 



xAikAH-'l^lf — (A.6) 

V U| k '- k ^ 6A{k 2 )6A(k 3 ) V ' 

It is worth mentioning here that the term proportional to A\ kept in |6) gives only an 0(L~ d ) contribution in the large 
box limit, which can be calculated exactly as 

9 ^f)V^|^^ (A.7) 

which is even smaller than the 0(L~ d/2 ) correction from J'x which we calculated before. 



Calculation of J 3: Substituting the definition (73 1 of a. we have 



Ji = Ewas (* + ^K**^* \^ws (^rvr^vr n k <% (a.8) 

£r(< 5 , 6>ffl)*L^sU=i + C2 «"» 3). (A.9) 
By symmetry, we only need to consider the first part, which has the following diagram pre-phase-averaged: 

34 



1 



> 



< 



2 



3' 



After phase averaging, there are three types of contributions. The type I diagrams are two leading contributions which 
contain two free wavenumbers. We show only one diagram here, the other obtained by exchanging 4 «-> 5 



The type II diagrams are four subleading contributions which contain one free wavenumber. Again we only show one 
diagram here, with mode 1 externally coupled: 



It turns out that the graphs with 1 externally coupled are the most important type II contributions to J3. The delta 
functions require lt4 = -ki and <x 4 = -1 in this diagram. Additional type II graphs can be obtained by symmetry. All 
other contributions are type III where all wavenumbers are pinned. 

Now consider large L. The amplitude factor changes the order of J3 by L~ 2d . In the term proportional 

to A\, the type I diagrams increase the order by I? d since they contain two free sums, giving an 0(1) contribution to 
STi- In the term proportional to the type I diagrams give zero contribution because they require fXj to vanish. The 
leading 0(1) contribution now comes from the type II diagrams with mode 1 externally coupled and fi\ = 1. The type 
III diagrams contain no free wavenumber and give contributions at most 0(Lr d ). The leading contribution to J3 in 
the large-box limit is thus: 





J3 




0-2,0-3 k 1,2,3 



0-2,0-3 1,2,3 k^kt.-ki 



xE T (a)i + cj-u coi - cr 2 co2 - c^ws^kLO^+rxjk, + ( 2 <^ 3) 



(A. 10) 



Substituting Jj 



(2nlL) d Jj and A\ = iA(k.{) and taking the limit L — > oo using (84i, one finds 




o;=(-l,o-2,o-3) 



35 



$ 2 Z 



<M(ki)<M(k 3 ) 

-9j]<WiVu [I <W Z ^ 2 f ^ 2 f d%8 d (g:-k)E T (co(ki) + coi-kila - co(^)) 
1 k#ki,-ki — -■> J J 



0"=(— l,Cr2,CT3) 



ffi>) + (2 - 3). 



(A. 11) 



Calculation of 3a,: Similarly, substituting ( 72 1 for a\ 



XA r (o-2<y2 + O-3W3 - <yi)A7-(cr 4 W 4 + cr 5 (L) 5 - W 1 )(5 kl>(J2k2+tr , k ,5k 1 , t r 4 k4+o-5k 5I 

which can be represented by the following diagram before phase-averaging: 

4\ A 



(A. 12) 




After phase averaging, there are two types of contributions. Type I diagrams have only the 1 mode externally coupled, 
corresponding to the following diagram and its partner with 2 <-> 3: 




These are formally the leading contributions, with one free wavenumber. However, the wavenumber delta functions 
require k; = and thus this graph vanishes identically. All other nonzero contributions are type II where there are no 
free wavenumbers. 

Now consider large L. The amplitude factor J\ ^[J^TyJ^Js in J~A changes the order by L~ 3d . The largest possible 
contribution comes from the term proportional to jjii^j - 1), which is a large factor of 0(L 2d ). However, such 

contributions are still only 0(L~ d ). Hence ?J\ is negligible in the large box limit. 

Calculation of J 5 : Write J 5 = £1*2 {/M2C81 + £2) + (M + j|)f (Si - S3)}, where J72I gives 



s 2 = 



k ' ^ 3,4,5,6 < k 

XA(cr 3 W 3 + O-4W4 - 0>i)A{o- 5 (l) 5 + CT 6 a»6 - W2)5 kl ^ 3k , +tr4k4 <5 k2jtr5k5+0 - 6 k 6 

n^w) = z c^^^kwwii'f 

k > $ 3,4,5,6 ' k 

XA(cr 3 W 3 + O-4W4 + dL» 1 )A(o-g ( y 5 + CT 6 a)(, - W2)^-kl,0-3k3+£r 4 k4^k2,0-5k5+Cr 6 k 6 
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' k 3,4,5,6 \ k 

XA(CT 3 W 3 + O-4W4 - Wi)A(cr 5 W 5 + cr 6 W6 + W2)<5k 1 , t r3k,+ t r4k4^-k 2 ,o- 5 k 5 + t r (1 k 6 - (A.13) 



Note under the interchanges 1 «-> 2, 3 <-> 5, 4 <-» 6, that S2 «-» S3. All three of the S's contain the same amplitude 
factor -\JJ1J2J3J4J5J6 that changes the order by L~ id . We list their diagrams as follows, 




After phase averaging, there are three types of diagrams for the S's. The type I diagrams are leading contributions 
with no external couplings. Shown here is one type I diagram for Si : 




with others obtained by 3 «-» 4 and 5 <-> 6. These contain two free wavenumbers. For S2, S3 there are additional type 

I contributions with 1 and 2 coupled, but these do not contribute to J's because of the restriction that The type 

II diagrams are subleading contributions which contain one free wavenumber. Shown here is one such contribution 
for Si with both 1 and 2 externally coupled: 




The delta functions here require k2 = -ki . In the corresponding diagrams for S2, S3 the delta functions instead require 
ki = k2, so these do not exist. (Contributions to S2, S3 with 1 and 2 coupled to the same external blob coincide with 
the type I diagrams that do not contribute to J5.) There are also type II contributions for the S's with only one of 1 
or 2 externally coupled. For example, one such contribution is obtained for T>\ by exchanging 1 «-» 4 in the above 
graph and then the delta functions require that k2 = <x 4 k4 so that k2 = — kt, <x 4 = -1. Similar type II contributions 
exist for S2, S3. Other type II contributions with neither 1 nor 2 externally coupled must have both /Ji = ^2 = and 
do not contribute at non-vanishing order to STs- Type III contributions with additional external couplings have no free 
wavenumbers. 

Now consider the size of the various terms for large L. For £1^2 ^1^2(^1 + S2), the leading contribution comes 
from type I diagrams and is 0(LT d ). For £ 1#2 ^(S, - S3), the type I diagrams give zero contribution since ji/2 = 
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and type II diagrams are at most 0(L d ). For 2i#2 &u type 1 diagrams give zero contribution and type II diagrams 
with both 1 and 2 externally coupled give an 0(1) contribution. The corresponding term 2 1^2 ajj'^ ~ 0(L^ d ) from 
type III diagrams, since, as noted above, there do not exist type II diagrams for S3 with 1 and 2 externally coupled to 
distinct blobs. In large box limit, the term STs is given, to leading order, by 

Js ~ [] ^» Z [ dl ' k2 f^V(£'k)A r (£' W (k))A T (^' W (k')) 

1 k#±k[ £=(- 1,0-2, 0-3 ) ^ ^ 



1 ' / 



where k = (ki,k2,k3) and k' = (— k1.k2.k3). Now making the substitution A^ = iXQs) and taking the limit, one 
obtains, 



(« ElWl J 5 ) ; ~ -^^,ii„i f] ^ Z fd d k 2 d%6 d (a • k)A r (<I • w(k))A r (-cr • w(k')) 

1 k#±ki £=(-,0-2,0-3) 



• (A. 15) 



Appendix B. Derivation of the PDF Hierarchy 

We here evaluate the JT-terms that contribute to A" (M) (/l,//, T). Most of the essential work has already been done 
in |Appendix A| As discussed in the text, the main difference is that mode 1 is now discrete for Ji - J4 and assumes 
only M values, while forJs modes 1,2 are both discrete. All other modes are continuous in the infinite-box limit. One 
must consider carefully whether free wavenumbers in graphical summations are discrete or continuous to see whether 
their contribution is O(M) or 0{L d ). 

\ The graphs contain no free wavenumbers, so that J7i is still 0(L~ d l 2 ). 

$2'. The leading order 0(1) contribution comes again from the type I diagrams, which contain two free wavenum- 
bers. One of these free wavenumbers is mode 1 which is discrete, but this is compensated by the extra factor of L d 
from the prefactors and A^J^, which now contribute equally. One thus obtains by substituting A^ = i(L/2n) d A m 
and J m = (2nlL) d J m the result 

j=\ £=(-1,0-2,0-3) 

x((iA j -A)j j )jy^ i ^) j 

M „ 

= -9i^,o J] J] J d%d%6 d ( Z -k j )\A T ( Z -cj(k))\ 2 \H0 

j=i oj=(-i,o- 2 ,o- 3 ) 

a a2 y(M+2) 

x(A j + A]i-) a _ _ . (B.l) 

J dAj 3A 2 dA 3 Uz=As=o 

Here k^. = (ky, k2,k3). 

J3: There are two contributions, from the two terms in the prefactor, Ak { and jUk,/27k, • For the A^ term, the 
leading 0(1) contribution comes from the type I diagrams, just as for ^2- One of the two free wavenumbers is discrete 
mode 1, but this is compensated by the 0(L d ) prefactor A^. For the fi^JU^ term, the leading 0(1) contribution 
comes from the type II diagrams with modes 1 and - 1 externally coupled, so that fi^ = ju_k, = 1 ■ The type II diagram 
has one free continuous wavenumber. Thus, its contribution is also 0(1) if the set of M wavenumbers includes at least 
one pair of opposite wavenumbers. Assuming this is not so, the final result is: 

M — , r 

(f ZlW 'j3)" ~ -18^,o 2 Z a2 J d%d%S?(a-^E T (0,gi-(o(^)\^ 

j-l (X-(-l,a"2,0"3) 
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x (iXjJjJ^e^^-) + (2 <-> 3). 
W)£ 2 ^2 f ^M%^(^-k.)£r(0,^-w(k.)) 

j=l a:=(-l,a-2,a- 3 ) 



M 

18i5, 



a2 -r(M+l) 

x^'lT/jrl— = _ + (2<->3). (B.2) 

1 rXlu'AI, U=o 



However, if the M wavenumbers include pairs of opposites, then there is an additional contribution 

M p p 

( e i*v* j3 y» „ -9jvVj, 11 5 a-o Z ^ \ dd ~ kl \ dd ~ k ^ d ^-h 

j=\ m*j-j gy=(- 1,0-2,0-3) ^ ^ 

T- ,0-2,0-3 rr+ ,0-2,0-3 / j;,U„l, M-l 



J J — -J k ; ,k 2 ,k 3 k„k 2 ,k 3 \ V J l j 

+(2<-»3). (B.3) 

J\. The type I contribution vanishes and the leading contribution from type II diagrams with no free sums remains 
0(L~ d ). 

J5: The amplitude factors contribute 0(L~ 3d ) while the prefactors are <9(L M ). Thus, to give an overall 0(1) 
contribution, the £ factors must contain a free continuous wavenumber. For Si the type I diagrams have two free 
discrete wavenumbers and type II diagrams with mode 2 only pinned have one free discrete wavenumber, but type 
II diagrams with both 1 and 2 pinned have one free continuous wavenumber. For S2,S'i, there are no graphs with 
modes 1 and 2 distinct, contributing to J5, that have a free continuous wavenumber. Type I diagrams have two free 
wavenumbers, both discrete if 1 and 2 are uncoupled. Type II diagrams with only one of the modes 1 or 2 externally 
coupled have one discrete free wavenumber. 

Thus, the only 0(1) contribution to J" 5 is from £1*2 {/tk,^k 2 + -7^ + B^lt}^ 1 anc * ^ diagrams for S\ 

with both 1 and 2 pinned. Corrections are at most 0(M/L d ). The 0(1) contribution is thus like that obtained before 
for J5, except that the mode-1 sum is over just M discrete values and all terms in the prefactor contribute equally. 
This gives 



7=1 mtl,-l ~1 1 ~ tr=(- ,1X2,0-3) 

xAAct ■ (o(k ))A T (-a ■ (o(K))H--™H + -™ I Jjjlj 2 he 1 ^ . (B.4) 

— —J — —J k,,k 2 ,k 3 -k,,k 2 ,k 3 \ > j 
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